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CHAPTER  I 


INTRODUCTION 


In  electromagnetic  problems  the  main  task  is  usually  to  deter¬ 
mine  the  current  distribution  on  the  surface  of  a  scatterer.  It  is 
a  relatively  simple  problem  to  obtain  other  parameters  we  desire  from 
a  knowledge  of  the  surface  current  distribution.  Obtaining  the  cur¬ 
rent  distribution  usually  means  we  must  solve  integral  equations  de¬ 
rived  from  Maxwell's  equations  and  the  boundary  conditions. 

One  of  the  well  known  methods  for  solving  integral  equations 

f 

is  the  moment  method  [6] ,  [12] ,  [16] ,  which  converts  the  integral  equa¬ 
tion  to  linear  algebraic  equations  having  as  many  unknowns  as  the 
number  of  basis  functions  used  to  approximate  the  surface  current 
on  the  scattering  body.  Even  though  the  moment  method  gives  us  rela¬ 
tively  accurate  results  for  arbitrarily  shaped  bodies,  its  practical 
use  is  usually  limited  to  bodies  which  are  not  large  in  terms  of  a 
wavelength. 

Another  powerful  method.,  which  is  known  as  the  geometrical 
theory  of  diffraction  (STD),  was  introduced  by  Keller  (1953).  GTD 
is  a  ray  optical  method  which  uses  the  leading  terms  of  the  asymp¬ 
totic  approximation  to  the  integral  of  interest.  Recently,  Kouyoumjian 
and  Pathak  [5]  have  developed  new  formulations  for  diffraction  coef¬ 
ficients  which  produce  the  better  results  in  the  transition  regions. 
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Usually  GTO  is  known  as  a  high  frequency  method  since  it  has  an  ad¬ 
vantage  when  applied  to  scatterers  which  are  electrically  large. 
However,  GTO  solutions  are  available  only  for  relatively  simple  can¬ 
onical  geometries. 

Although  both  methods  considered  above  provide  us  with  relatively 
accurate  results  in  their  own  applications,  there  are  many  electro¬ 
magnetic  problems  which  cannot  be  solved  by  either  method  alone.  To 
solve  such  problems  by  using  to  advantage  the  two  methods  mentioned 
above,  recently,  hybrid  techniques  which  combine  them  have  been  de¬ 
veloped.  One  approach,  which  has  been  introduced  by  Thiele  and 
Newhouse  [2]  in  1975,  combines  moment  methods  with  the  GTD  ( i . e. , 
MM-GTD  techniques).  This  technique  extends  the  moment  method  through 
the  use  of  the  GTD  to  obtain  a  new  modified  impedance  matrix  for  the 
moment  method  solution.  The  application  of  this  method  was  demon¬ 
strated  by  solving  the  problems  such  as  a  monopole  near  a  conduct¬ 
ing  wedge,  a  monopole  at  the  center  of  a  flat  plane  or  a  circular 
disc,  and  a  monopole  near  a  conducting  step.  A  little  later  in  the 
same  year,  Burnside  and  his  colleagues  [3]  showed  another  hybrid 
approach  by  combining  the  GTD  with  the  moment  method  ( i . e. ,  GTD-MM 
technique).  In  this  technique,  the  diffraction  coefficient  is  treated 
as  an  unknown  coefficient  in  the  moment  method  and  was  determined 
by  numerical  techniques.  The  application  of  the  method  was  demon¬ 
strated  by  obtaining  the  surface  currents  on  a  perfectly  conducting 
wedge,  and  square  and  circular  cylinders.  Later  Mittra  [7], 18],  also 
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introduced  a  new  approach  for  combining  the  integral  equation  and 
the  GTO  by  using  the  fact  that  the  Fourier  transform  of  the  unknown 
surface  current  is  proportional  to  the  scattered  far  field  and  dem¬ 
onstrated  its  applications  to  a  perfectly  conducting  infinite  strip 
of  zero  thickness,  a  thin  conducting  rectangular  plate,  and  later 
to  a  circular  cylinder  in  two  dimensions. 

The  GTD-MM  method  which  combines  the  GTD  with  the  moment  method, 
requires  a  priori  knowledge  of  the  asymptotic  form  of  the  current 
away  from  the  diffraction  or  moment  method  region.  Therefore,  the 
GTD-MM  technique  considered  above  is  difficult  to  apply  to  arbitrar¬ 
ily  shaped  bodies  where  the  GTD  current  forms  are  not  known  a  priori. 

In  this  paper  a  new  hybrid  method  which  overcomes  the  disadvan¬ 
tage  mentioned  above  is  presented.  This  method  which  combines  the 
moment  method  with  an  asymptotic  technique  as  orginally  introduced 
by  Molinet  Cl] ,  does  not  need  any  priori  knowledge  of  current  forms 
away  from  the  moment  method  regions.  Thus,  the  hybrid  method  of  this 
paper  offers  an  advantage  in  solving  for  currents  in  GTD  transition 
regions  where  the  form  of  the  current  is  often  difficult  to  determine. 
Further,  this  hybrid  method,  which  uses  the  geometrical  optics  current 
as  the  dominant  contributor  to  the  initial  current,  works  surprisingly 
well  even  for  extremely  small  bodies  in  the  Rayleigh  region.  One 
would  expect  this  technique  to  work  well  for  large  bodies  since  it 
uses  the  geometrical  optics  current  which  is  asymptotic  with  respect 
to  frequency.  The  fact  that  it  also  works  well  for  low  frequency 
is  indeed  an  unexpected  bonus. 
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The  integral  equations  originally  introduced  for  this  new  hybrid 


method,  were  derived  for  the  wedge  problem  and  are  not  suitable  for 
other  shaped  bodies.  Therefore,  in  Chapter  II,  we  consider  the  gen¬ 
eral  theory  for  the  method  and  also  derive  the  general  integral  equa¬ 
tions  which  are  suitable  for  arbitrarily  shaped  bodies.  In  the  de¬ 
rivation  of  the  equations  the  asymptotic  currents  are  assumed  to 
approximately  represent  the  dominant  currents  on  the  surfaces  of  a 
conducting  body.  Then  the  difference  current  between  the  approximate 
and  exact  currents  is  obtained  by  using  the  moment  method  in  the 
moment  method  regions  which  are  formed  around  the  shadow  boundary 
regions  or  around  sharp  geometrical  discontinuities.  In  the  next 
step,  we  determine  the  currents  in  the  non-moment  method  regions  by 
obtaining  the  contribution  to  the  surface  currents  induced  by  the 
currents  in  the  moment  method  regions.  The  process  may  be  repeated 
to  obtain  higher  order  solutions. 

Applications  of  the  method  to  a  perfectly  conducting  wedge,  square 
and  circular  cylinders  in  two  dimensions  are  considered  in  Chapters 
III,  IV,  and  V,  respectively.  In  Chapter  VI,  extension  of  our  hy¬ 
brid  method  to  the  three  dimensional  geometry  is  also  considered  by 
using  a  perfectly  conducting  sphere  as  an  example.  The  results  of 
the  surface  currents  on  a  wedge  and  a  square  cylinder  were  compared 
with  the  solutions  from  the  GTO  whose  diffraction  coefficients  are 
based  on  the  results  derived  by  Kouyoumjian  and  Pathak  [5].  On  the 
other  hand,  for  a  circular  cylinder  and  a  sphere,  the  comparisons 
were  made  with  the  exact  eigenfunction  solutions. 

f 
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Usually  for  bodies  having  only  plane  surfaces,  the  general  in¬ 
tegral  equations  for  the  method  reduce  to  relatively  simple  forms, 
but  not  for  bodies  composed  of  curved  surfaces.  Also  observing  those 
equations  for  the  plane  surfaces,  we  can  realize  that  singularity  prob- 

t 

lems  are  not  involved  in  them.  However,  this  problem  will  arise  for 
curved  surfaces  and  must  be  considered  for  more  accurate  results. 
Theoretically  when  we  evaluate  the  principle  value  of  integration, 
we  must  exclude  the  singular  point  for  numerical  integration,  but 
practically  we  exclude  the  patch  of  a  finite  area  containing  the  sin¬ 
gular  point  instead  of  just  the  singular  point.  We  see  that  this 
practical  calculation  will  increase  errors  as  long  as  the  surface 
is  not  divided  into  very  many  subpatches  so  that  the  patch  areas  be¬ 
come  very  small.  However,  considering  the  computer  running  time  and 
memory  storage,  it  may  not  be  efficient  to  divide  the  surface  into 
so  many  divisions.  Therefore,  in  Chapters  V  and  VI  we  derive  some 
constants  which  will  roughly  provide  the  contribution  lost  by  exclud¬ 
ing  the  finite  patch  area  instead  of  a  single  point  for  the  evaluation 
of  the  principle  value  of  integration,  and  the  effect  of  those  con¬ 
stants  is  shown  by  comparing  the  results  with  and  without  them.  By 
using  those  constants  in  the  integral  equations  for  the  method,  we 
can  decrease  the  number  of  divisions  for  numerical  integrations  and 
thus  computing  time  to  obtain  the  same  results  that  can  be  obtained 
by  increased  number  of  divisions. 
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In  the  examples  which  are  considered  in  Chapters  III,  IV,  V, 
and  VI,  the  impedance  matrix  for  the  moment  method  is  derived  exactly 
for  the  wedge,  but  approximately  for  the  square  and  circular  cylinders, 
and  the  sphere.  Even  though  we  can  derive  the  more  accurate  expres¬ 
sion  for  it,  it  may  not  be  useful  for  numerical  calculations  because 
of  its  complexity. 

In  Chapter  VII,  summary  and  conclusions  are  presented  along  with 
some  of  the  advantages  and  disadvantages  of  this  method. 
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CHAPTER  II 


GENERAL  THEORY 


2.1  First  Order  Approximation 

Let's  consider  an  arbitrary  shaped  3-dimensional  perfectly 
conducting  body  illuminated  by  a  plane  wave  as  shown  in  Figure  2.1. 

We  divide  the  total  surface  S  of  the  body  into  and  Ss  which  represent 
the  lit  and  shadowed  regions,  respectively.  The  separation  boundary 
of  S  is  the  shadow  boundary  which  is  shown  by  the  line  AB  in  Figure 
2.1.  In  addition,  we  split  each  of  these  regions  into  two  regions, 
a  moment  method  region  (MM-region)  and  an  asymptotic  region  (AS-region) 
which  are  represented,  respectively  by  S^M  and  S^s  in  the  lit  region 
or  by  and  S^s  in  the  shadow  region.  Therefore,  we  can  write  the 
following  expressions  regarding  the  surfaces. 

?  ■  *  is  <2J> 

sS  =  SMM  +  SAS  •  (2*2) 

Usually  the  moment  method  region  is  smaller  than  the  asymptotic  region. 

Next  let  the  total  current  be  denoted  by^'(^)  in  the  lit  region 
and  by  in  the  shadowed  region.  The  current  can  then  be  expressed 

in  terms  of  JM^j{R)  and  3*^(1?)  or  3^(R)  and  ^$(R)  as 
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,£  /  /  !  /  I 

‘AS  /#  /  //  /  / 

/  Ws  //  /  / 

/  y/  SMM/  7  ' 

_  /  , r  -V'  /  / S8 

/  /  /  °as 


U-'-L'-p-r 


z  Figure  2.1.  Three  dimensional  surface  is  divided 

into  four  regions. 

=  jLflO  +  7*1(10  (2.3) 


?S(ff)  =  ^(R)  +  ^(*)  (2.4) 

where  U^ff)  and  jj^ff)  are  the  total  currents  on  S^M  and  S^M,  respect¬ 
ively,  and  3^(R)  and  ”3^(R)  represent  the  total  currents  on  S^s  and 
S^,  respectively. 

We  will  use  the  usual  form  of  Maue's  magnetic  field  integral 
equation  to  derive  the  general  equations  for  our  hybrid  method.  The 
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current  at  an  observation  point  P  (see  Figure  2.1)  will  satisfy  the 
magnetic  field  integral  equation  given  by 

3(1?)  =  2nxHi  (1?)  +  2nx/  J(TT* )  x  V'G(r)ds'  (2.5) 

s 

where  1?  and  1?'  are  the  position  vectors  of  the  observation  point  P 
and  the  source  point  P',  respectively,  and  r  *  ll?-!?' |  is  the 
distance  between  the  points  P  and  P*.  The  unit  vector  n  is  outward 

and  normal  to  the  surface  at  P.  H*(J?)  is  the  incident  magnetic  field 
at  P  on  the  surface.  The  symbol  f  represents  the  principle  value  of 
integration. 

The  free  space  Green's  function  G(r)  for  the  two  dimensional 
problem  is  given  by 

G(r)  *  h  H£2)(Br)  (2.6) 

and  for  the  three-demensional  problem  is  given  by 
“*  jB  r 

( 2 . 7 ) 

where  H^(Br)  is  the  zero  order  second  kind  Hankel  function,  and  $  is 
the  free  space  propagation  constant.  Now  v'G(r),  where  the  prime 
refers  to  the  coordinate  system  of  the  source  points,  can  be  obtained 
as 

7'G(r)  =  -  Hj2)(6r)  r'  (2.8) 
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for  the  two-dimensional  problem  and 


_  j  py* 

V’G(r)  =  -  (j8  +  -p)  r'  (2.9) 

for  the  three-dimensional  problem,  where  ?'  is  the  unit  vector  in  the 
direction  of  IT'-IT,  and  H^(gr)  is  the  first  order  second  kind  Hankel 
function.  The  time  dependent  factor  eJU)t  is  suppressed  throughout. 

The  integration  surface  S  in  Equation  (2.5)  can  be  divided  into 
£  and  Ss,  and  according  to  those  integration  regions,  the  integral 
can  be  separated  into  two  parts.  For  observation  points  in  the  lit 
region,  the  current  ^(R)  will  satisfy  Equation  (2.5)  or  the  following 
equation. 

=  2nxHi(ff)  +  Zf\  x  ffti R')  x  7'G(r)ds' 

S* 

+  2ft  x  /  JS(R')  x  v'G(r)ds*  (2.10) 

Ss 

The  second  and  third  terms  are  the  currents  induced  by  the  currents 
at  source  points  in  the  lit  and  shadowed  regions,  respectively. 

Since  contribution  of  the  current  3S(R')  is  usually  small  com- 

O  _ 

pared  to  that  of  "J  (R ' ) ,  the  third  term  of  Equation  (2.10)  is  small 
compared  to  the  second  term,  and  also  both  of  them  are  small  compared 
to  the  physical  optics  term  2nxH1(ff).  Therefore,  ignoring  the  third 
term  in  Equation  (2.10),  we  can  make  the  approximation  of 
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(2.11) 


3*flD  *  3jp{R) 

where  the  optics  region  current  Jgp(R)  is  given  by 

3j  00  =  2nxfi1  (I)  +  2nx  f  (**)  x  7’G(r)ds'  (2.12) 

S®' 

Also  we  can  write  the  following  integral  equation  for  observation 
points  on  the  shadowed  region. 

1S<$)  =  2nxHi(ff)  +  2nx  f  ^(K*)  x  v’G(r)ds' 

S* 

+  2nx  f  JS(R ' )  x  VG(r)ds '  .  (2.13) 

Ss 

If  we  replace  (R 1 )  in  Equation  (2.13)  by  the  approximate  current 

n  _ 

1  m.  which  is  given  in  Equation  (2.12),  then  the  current  "3s (1?)  is 
no  longer  exact,  but  instead  is  approximate.  If  we  represent  this 
approximation  by  ^p(^),  then  Equation  (2.13)  becomes 

3jp(10  -  2nxRi (IT)  +  2nx  /  ^p(R‘)  x  v'G(r)ds' 

Sl 

+  2nx  /  ^p(R')  x  V'G(r)ds'  .  (2.14) 

Ss 

Since  JqP(I)  in  Equation  (2.12)  and  J^p(R)  in  Equation  (2.14) 
are  approximations  for  J^R)  and  J^R),  respectively,  additional  cur¬ 
rents  are  required  to  determine  the  correct  currents  on  both  surfaces. 
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Let  these  additional  currents  be  I* (ft)  and  Ts(ft)  on  the  surfaces  S8, 
and  Ss,  respectively.  Then  the  exact  currents  in  both  the  regions 


*  (R>  *  3*  <*)  +  T^TT) 


(2.15) 


3*00  -  3^00  *  TS(R) 


(2.16) 


According  to  Equations  (2.3)  and  (2.4),  the  currents  T^ft)  and  Ts(ft) 
can  be  expressed  as 


Tl(R)  =  T‘ <R)  ♦  T*  (R) 


(2.17) 


IS(R)  •  IjL(R)  ♦  lL(R) 


(2.18) 


Then  from  Equations  (2.3),  (2.4),  (2.15),  (2.16),  (2.17),  and  (2.18) 
we  can  express  the  total  currents  on  the  regions  S^M,  S^M,  S^,  and 

SAS  aS 


J 1(R)  ■  3*  (R)  +  li(R) 


3L(R)  -  3V(R)  *  I*  Of) 


3f.(R)  ■  3*  (R)  ♦  T*  (R) 


(2.19) 


(2.20) 


(2.21) 
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^AS<*>  *  V>  *  ^S<*>  • 


(2.22) 


The  currents  at  observation  points  on  the  moment  method  region 
also  will  satisfy  Equation  (2.5).  Hence, 

*  ZffxJT'flT)  +  2nx  f  x  V6(r)ds' 

Sl 

+  2nx  /  ^(ff')  x  V'G(r)ds'  (2.23) 

Ss 

and 

4flf)  =  2nxH1(E)  +  2nx  /  ^(fT1)  x  V'G(r)ds' 

Sl 

+  2nx  f  3S(R‘‘)  x  V'G(r)ds'  .  (2.24) 

Ss 

Substituting  Equations  (2.15),  (2.16),  (2.19),  and  (2.20)  into 
Equations  (2.23)  and  (2.24),  and  then  using  Equations  (2.17)  and  (2.18), 
we  obtain  the  following  two  equations. 

tyR)  +  14,00  =  2nxHi  (IT)  +  2nx  f  3jp(ff')  x  V'G(r)ds' 

S* 
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+  2nx  /  U*  (ff*)  x  V'GCrJds1  +  2nx  /  T*  (*')  x  V'G(r)ds' 


+  2nx  /  Tfe(E')  x  V'G(r)ds' 


(2.25) 


ZDi*)  +  ijLflD  =  2fixH1(ff)  +  2nx  f  J*(ff')  x  VG(r)ds 


+  2nx  j  Tmm(^')  x  V’G(r)ds'  +  2nx  f  iL(E')  x  V*G(r)ds * 


+  2itx  /  3*  (R* )  x  V'G(r)ds'  +  2nx  f  TL(R')  x  V*G(r)ds' 


+  2riV  /  ljjs(ff')  x  y*G(r)ds'  . 


(2.26) 


Subtracting  Equation  (2.12)  from  Equation  (2.25)  and  Equation 
(2.14)  from  Equation  (2.26),  we  obtain  the  currents  in  the  lit  and 
shadowed  moment  method  regions  which  are  given  by 


‘‘**',**-v 


=  2nx  /  TLtf^xV'GfrJds'  +  2nx  f  TjLOT' )xV'G(r)ds' 


+  2fix  /  r _(R' )xy G(r)ds*  +  2fix  /  I~<.(R' )xV'G(r)ds' 


+  2ftx  f  TjLfR'Jxv'GfrJds' 


(2.27) 


TL(I)  «  2fix  /  IuM(R' )xV'G(r)ds'  +  2nx  /  I*  (If' )xV'G(r)ds' 


+  2 fix  f  T^^'Jxv'GCrJds'  +  2nx  f  )x7’G(r)ds 1 


(2.28) 


Since  there  are  four  unknown  currents  TiW*>-  4<*>.  TAlsro, 
and  lj^(R:)  in  two  equations,  we  cannot  solve  for  them  without  making 
proper  approximations.  In  Equation  (2.27)  it  is  clearly  seen  that 
the  current  term  containing  ^(R1)  Is  dominant.  However,  it  is  hard 
to  compare  the  other  current  terms  with  each  other  analytically.  As 
we  know,  the  currents  on  the  right  side  of  Equation  (2.27)  represent 
the  currents  Induced  by  3*p(R)  and  the  currents  in  the  moment  and 
asymptotic  method  regions. 


Now  considering  that  observation  points  are  limited  to  the  moment 
method  region  and  that  the  current  in  the  nearest  region  to  the  obser¬ 


vation  points  will  most  contribute  to  the  current  at  the  observation 
points,  we  can  ignore  those  current  terms  containing  or 

in  both  Equations  (2.27)  and  (2.28)  to  derive  the  first  order  approxi¬ 
mations  for  T|4,(R)  and  Tj^ff).  Therefore,  the  first  order  approxi¬ 
mations  for  T^(E)  and  T^(W)  can  be  written  as 

I|4,(R)  =  Znx  /  Tjjjn,(R' )xV'G(r)ds'  +  2nx  /  TjM(R' )xV'G(r)ds' 
SjL  Sum 

MM  MM 

+  2nx  f  Jq  (R'JxV'G^ds'  (2.29) 

_  c 


and 


Tj^ff)  =  2nx/  )xV'G(r)ds'  +  2nx  f  1^,(1? ' )xV'G(r)ds* 

si,  sL 

MM  MM 

(2.30) 

We  should  realize  that  without  ignoring  those  current  terms  containing 
Tj-flP )  or  T^j(R')  in  Equations  (2.27)  and  (2.28),  we  cannot  derive 
equation  sets  which  can  be  solved  by  the  moment  method.  On  the  other 
hand,  since  the  current  ^(E-')  is  known,  we  can  solve  Equations 
(2.29)  and  (2.30)  simultaneously  using  the  moment  method.  By  the  same 
procedure,  we  can  derive  integral  equations  for  IA5(ff)  andTAS(ff). 
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However,  we  can  write  down  directly  those  equations  by  noting  that 
the  currents  on  the  asymptotic  regions  also  satisfy  Equations  (2.29) 
and  (2.30)  for  their  first  order  approximations.  Thus,  we  can  write 
the  following  equations  for  T^(ft)  and  '«<*>• 

1 

T}s(ff)  -  2nx  /  l^')x7'G(r)ds*  +  2fix  f  !*,(*•  )xv*G(r)ds' 

SMM  SMM 

+  2nx  f  3!  (JT'Jxv'GtrJds'  (2.31) 

Ss 

and 


(R)  =  2flx  j 
S®’ 

Juu 

nn 


JMM(R ' )xV*G(r)ds  *  +  2nx 


/  iLOl'JxV'GfrJds1 

Ss 

5mm 


(2.32) 


T^<-(R)  and  Tj^(ff)  can  be  obtained  simply  by  substituting  the 
pre- determined  currents  4(1*)  and  TJmOT)  into  Equations  (2.31)  and 
(2.32). 

The  above  first  order  approximations  may  not  be  enough  to  give 
the  correct  surface  currents  on  some  complex  bodies.  Therefore,  we 
will  derive  equations  for  the  second  order  approximations  for  TjJ^(R) 
and  Tj^|(ft)  in  the  next  section. 
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Second  Order  Approximation 


We  can  obtain  the  second  order  approximations  for  T^M(E’)  and 
by  substituting  Equations  (2.31)  and  (2.32)  into  Equations 
(2.27)  and  (2.28),  which  are  given  by 

TJjmO?)  *  ^(ff)  +  2n£x  /  [  2n£x  f  T^M(^")xv,,G(r1)ds"] xV'G(r2)ds* 

^AS  ^MM 

+  2n£x  f  (2n£x  f  TjM(R")xV"G(r1)ds" ]xv'(r2)ds' 

pJl  cS 

"'AS  ^MM 

+  2n£x  /  (2n$x  f  TjM(R")xv,,G(r1)ds" ]xV*G(r2)ds' 

pS  pj, 

5as  5mm 

+  2n£xf  C2ft$x  f  T^M(R")xv,,G(r1)ds"]x7,G(r2)ds' 

pS  pS 

^AS  *MM 

+  2rT£x  $  [2fi£x  /  j^p(lT'1)xvMG(r1)ds't]X7,G(r2)ds'  (2.33) 
S*AS  S$ 
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TjfjW  aTl^}  +  2ftsx  /  £2V  *  Tj5jjl(|(fi',)x7HG(rj)dsM]x7'G(r2)ds' 


+  2fisx  f  [2fiax  f  lJM(^M)x7"G(r1)ds,,]x7,G(r2)ds' 


+  2fi$x  f  [2n$x  f  I§M(^,')xv,'G(r1)ds,,]x7'G(r2)ds' 


+  2n$x  f  t2n$x  f  TJM(R,,)x7*,G(r1)ds,•]x7,G(r2)ds, 


Ss 

«S  ^MM 


+  2fi$x  f  [2nax  /  -fQ  (fM)x7"G(r1)ds"]x7,G(r2)ds'  (2.34) 


^AS  5 


where  I^(ff)  and  T^(£)  are  the  first  order  approximations  for  T^(JT) 

and  Tj^ff),  respectively,  and  n  and  ns  are  the  unit  normal  vectors 
in  the  lit  and  shadowed  surfaces,  respectively.  Also  r1  and  r2  are 
given  by 


r1  = 

r2  =  |R-f  | 


(2.35) 

(2.36) 


nmm 


The  second  order  approximations  for  T4s(ff)  and  T^s(tf)  have  the 
same  expressions  as  those  for  Tk<R> and  ■&<*>  ,  respectively,  except 
that  observation  points  are  moved  from  the  moment  method  regions  to 
the  asymptotic  regions. 

If  we  repeat  the  same  procedure,  we  can  obtain  equations  which 
will  give  us  higher  order  approximations  for  Tj^lf)  and  *  ^u" 

merically  this  procedure  can  be  achieved  easily  by  iterative  methods. 
2.3  Summary  and  Discussion 

We  have  derived  all  equations  for  the  first  and  second  order 
approximations  that  are  necessary  for  the  method.  For  summary  and 
later  convenience,  we  will  write  those  equations  for  the  first  order 
approximations  here  again 

■jJp(R)  =  2nxHi (IT)  +  2nx  f  )xv'G(r)ds'  (2.37) 

Sl 

^p(ff)  «  2nxHi(ff)  +  2nx  f  )xv‘6(r)ds' 

S4 

+  2nx  i  j*  (ff'hV’GfrJds'  (2.38) 

c 

Ss 
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^(10  =  2nx  /  )xV'G(r)ds '  +  2fix  /  f^R'  )xV'G(r)ds' 

et  cS 

mm  bm 

+  2nx  f  J^OT'Jxv'GfrJds'  (2.39) 

Ss 

2nx  /  ' )xV'G(r)ds'  +  2nx  /  T^ff' )xv'G(r)ds* 

C®’  cs 

nM  bMM 

(2.40) 

The  Green's  function  has  an  r~^  singularity  as  the  observation 
point  P  approaches  the  source  point  P'  on  the  surface  of  the  body. 

Since  the  second  term  of  Equation  (2.37)  and  the  third  term  of  Equation 
(2.38)  become  zero  if  P  and  P*  are  on  the  same  planar  surface,  the 
singularity  problem  is  not  involved  in  the  planar  surface.  However, 
there  will  be  as  many  singular  points  as  observation  points  for  a 
curved  surface  such  as  in  a  circular  cylinder  or  a  sphere. 

Theoretically  the  principle  value  of  integration  is  evaluated 
by  excluding  the  singular  point,  but  practically  it  is  done  by  exclud¬ 
ing  the  finite  patch  area  containing  the  singular  point.  Therefore, 
this  exclusion  of  the  patch  may  cause  some  errors  in  the  surface  cur¬ 
rent  calculations.  This  problem  will  be  considered  in  more  detail 
in  Chapters  V  and  VI. 
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fl  e 

When  we  use  the  iterative  method  to  solve  for  U*p(ff)  and  3op(E) 
in  Equations  (2.37)  and  (2.38),  we  have  to  consider  convergence  of 
the  result.  We  cannot  prove  mathematically  that  the  iterative  method 
produces  a  convergent  result,  but  instead  we  must  insure  it  by  numeri¬ 
cal  calculation. 

Also  it  appears  that  there  are  no  distinct  ways  to  determine 

% 

the  extent  of  the  moment  method  region  in  practical  numerical  calcu¬ 
lations.  However,  by  numerical  results  it  will  be  seen  that  the  extent 
of  the  moment  method  region  does  not  affect  the  surface  current  sig¬ 
nificantly.  Therefore  we  can  use  an  arbitrary  size  for  the  MM-region 
according  to  the  size  of  the  body.  This  is  very  desirable  for  prac¬ 
tical  applications  of  this  method,  because  if  the  currents  are  much 
dependent  on  the  extent  of  the  MM-region,  the  method  will  not  be  reli¬ 
able.  Usually  for  a  finite  dimensioned  body,  from  one  third  to  half 
of  the  size  of  the  body  will  be  reasonable  for  the  moment  method  re¬ 
gion.  If  we  increase  the  MM-region,  we  may  get  better  results,  but 
we  have  to  sacrifice  the  computer  running  time.  Also  there  will  be 
a  limit  to  increasing  the  MM-region  by  the  limit  of  computer  memory 
storages. 


22 


CHAPTER  III 


PERFECTLY  CONDUCTING  WEDGE 
3.1  Integral  Equations 

Consider  the  wedge  geometry  In  Figure  3.1,  where  the  one  side 
of  the  wedge  is  illuminated  by  a  TE  plane  wave  and  the  other  side  is 
shadowed.  We  set  up  the  moment  method  regions  around  the  edge,  which 
are  represented  by  C^  and  Cj^|  in  the  lit  and  shadowed  regions,  re¬ 
spectively.  For  the  wedge  problem,  the  currents  ^p(lT)  and  "3^p(lT) 
in  Equations  (2.12)  and  (2.14)  can  be  reduced  to 

7^(lf)  *  2nxH1(K)  (3.1) 

and 

3jp(*)  =  2nxfli(K)  +  2nx  /  ^p(K'  )xVG(r)dl'  .  (3.2) 

Cl 

Also  the  first  order  approximate  currents  in  the  MM-region  can  be 
simplified  as 


f 


and 


TMm(*)  =  2"x  f  Tj^(tf'  )x7'G(r)d£'  .  (3.4) 

r& 

UMM 

The  expressions  for  iJs^)  and  are  exactly  same  as  Equa¬ 

tions  (3.3)  and  (3.4),  respectively,  except  that  the  observation  points 
are  moved  to  the  asymptotic  region.  Therefore,  it  is  not  necessary 
to  write  those  equations  here. 

Now  let  us  consider  the  case  where  both  sides  of  the  wedge  are 
illuminated  by  a  plane  wave  as  in  Figure  3.2.  For  this  geometry,  the 

currents  (IT)  and  ^L(ff),  where  A  and  B  stand  for  the  two  surfaces 
op  op 

of  the  wedge,  can  be  expressed  as 

^p(ff)  *  2nxH1(tf)  +  2nx  /  ^p(^’  )xV,G(r)dJt*  (3.5) 

C8 

3^p(ff)  =  2nxH1(tf)  +  2nx  j  jjp(ff’ )x7'G(r)dr  .  (3.6) 

CA 

The  exact  current  on  the  surface  A,  can  be  obtained  from 

"3^x(1f)  *  2fixH*  +  2nx  /  ^x(]F')x7'G(r)di' 

CB 

+  2nx  /  ^x(^)x7’G(r)djl'  .  (3.7) 
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Figure  3.1.  Wedge  illuminated  by  a  TE  plane  wave,  where  ...one  side 
is  lit  and  the  other  is  shadowed. 

t 


Figure  3.2.  Wedge  illuminated  by  a  TE  plane  wave,  where 
both  sides  are  lit. 


However,  since  the  last  term  of  Equation  (3.7)  is  identically  zero. 
Equations  (3.5)  and  (3.7)  are  identical  except  for  the  subscript  on 

U(?[).  Thus,  the  currents  (If)  and  ^(ff)  represent  the  total  exact 
currents,  respectively  on  the  surfaces  A  and  B.  By  the  same  argument, 
this  is  also  true  for  Equation  (3.6).  This  means  that  we  don't  need 
any  additional  currents  (i.e.,  the  moment  method  region  currents). 
Actually  for  this  case  the  hybrid  method  is  not  used,  but  instead  only 

the  iterative  method  will  be  used  to  obtain  ^  (ID  and  J8  (1?). 

op'  '  op'  ' 

3.2  Impedance  Matrix 

To  obtain  the  impedance  matrix  for  the  moment  method  region  cur¬ 
rent,  we  will  express  Equations  (3.3)  and  (3.4)  in  scalar  forms  as 


(3.8) 

and 

&<*>  *2f  JA|(  >Gw^  )djt*  (3.9) 

C* 

m 

where 


GwU,f)  =  -  ft  Hj2)(Br) 


r 


(3.10) 
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and 


r 


J 


l2  +  l'2  -  2«' 


cos^ 


(3.11) 


The  above  result  can  be  readily  obtained  from  the  geometry  in  Figure 
3.3.  Substitution  of  Equation  (3.9)  in  Equation  (3.8)  leads  to 

I^U)  =  4  f  f  lMMU,)Gw(4%£,)Gw0l,r)dji'dr 
cs  r!> 

VIUJ  Luu 

rwi  nn 

+  2f  0®pU')GwU,i*)di'  .  (3.12) 

Cs 


Figure  3.3.  Geometry  showing  relationship  between 
su  and  $w  in  a  wedge. 
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k 


X  ‘><W  ■ 4 1  <*■) 

n- 1  pS 

UMM 

=  2/  *  (3.16) 

Cs 

Expressing  Equation  (3.16)  in  matrix  form 
[Zjn][lj]  =  [VjJ],  m,n«l,2 . N.  (3.17) 

allow  us  to  identify  the  elements  of  the  Impedance  matrix  [Z^l  as 

zmn  '  P<W  -  4  /  <3-'e> 

cs 

lmm 

where  m,n*l,2,3,...,N 

Usually  the  second  term  of  Equation  (3.18)  will  be  small  compared  to 
the  first  term.  The  elements  of  the  voltage  matrix  [V^l  are  determined 
by 

Vm  =  2  /  JopU')Gw(  Vt,d1,  ’  ( 3 . "1 9 ) 

Cs 
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3.3  Numerical  Results 


We  found  some  limitations  in  application  of  this  method  to  a 
wedge  after  having  tried  several  wedges  with  the  various  wedge  angles 
and  wave  incident  directions.  It  will  be  seen  in  the  examples  which 
follow  that  if  the  wedge  angle  (<b  )  or  the  wave  incident  angle  (*.) 

W  I 

is  small,  a  large  integration  limit  (XM)  is  required  to  obtain 
acceptable  results.  For  example,  if  <j»  is  90°  and  <j>.  is  30°,  then 

W  i 

XM  =  60A  is  needed  for  a  desirable  result.  Howevever,  if  we  decrease 
^  to  10°,  we  cannot  expect  a  desirable  result  even  with  a  value  of 
XM  larger  than  100a.  This  means  that  for  the  grazing  incident  wave, 
it  is  not  possible  to  obtain  an  acceptable  result.  Also  we  see  that 
the  wedge  with  very  small  wedge  angle  meets  the  same  difficulty. 

These  limitations  do  not  limit  the  method  theoretically,  but 
do  impose  a  practical  limitation  for  the  wedge  geometry.  Since  we 
deal  with  finite  structures  in  actual  problems,  these  limitations  may 
not  appear  in  other  practical  problems. 

We  obtained  the  surface  currents  on  the  lit  and  shadowed  regions 
up  to  5.0x  from  the  edge.  The  5.0A  distance  was  divided  into  100  sub¬ 
intervals  for  numerical  integration,  and  10  pulse  basis  functions  were 
used  in  the  0.5 A  wide  moment  method  region.  All  the  results  were  com¬ 
pared  with  independent  GT0  solutions.  Figure  3.4  shows  the  results 
for  4>w=<fcj=90°.  It  is  seen  that  the  magnitude  and  phase  of  the  surface 
currents  on  both  regions  are  in  very  good  agreement  with  the  GTD 
solutions.  In  Figure  3.5,  which  is  for  =90°  and  <J>,=30°,  we  can  recog- 
nize  a  little  difference  in  the  magnitude  of  the  shadow  region  current 
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in  the  region  from  4. Ox  to  5. Ox.  Note  that  this  result  was  obtained 
with  XM=60X  and  about  thirty  minutes  of  computing  time.  The  results 
for  the  surface  currents  for  the  geometry  in  which  both  surfaces  are 
illuminated  by  the  incident  wave  are  shown  in  Figure  3.6.  In  this 
case,  as  we  mentioned  before,  the  moment  method  region  current  was 
not  used  to  obtain  the  total  surface  current.  The  agreement  with  the 
GTD  solutions  are  shown  to  be  very  good  except  in  the  region  less  than 
0.2X.  We  show  the  results  for  *=120°  and  <j> .  =30°  in  Figure  3.7  and 
quite  good  agreements  with  the  GTD  results  are  indicated  except  near 
5. OX.  However,  for  the  reverse  geometry  (i.e.,  for  <t>w=30°  and  ^  =  120°) 
the  results  become  worse  as  shown  in  Figures  3.8  and  3.9.  The  results 
in  Figures  3.8  and  3.9  were  obtained  with  the  extent  of  the  moment 
method  region  (CMM)  of  0.5X  and  1.0X,  respectively.  As  was  expected, 
better  results  were  obtained  with  the  larger  CMM.  Even  though  many 
discrepancy  points  are  observed  on  the  magnitude  curve  for  0^=0. 5X 
the  results  are  still  acceptable  except  in  the  region  very  near  the 
edge.  Of  course,  better  results  can  be  obtained  with  increased  XM 
in  both  the  geometries  considered  above. 

By  observing  Figures  3.10  and  3.11  which  are  for  CMM=0.25X  and 
1 .Ox,  respectively,  we  can  realize  the  effect  of  the  size  of  CMM. 

Those  results  show  that  the  extent  of  the  moment  method  does  not  affect 
the  surface  currents  significantly.  As  we  increased  or  decreased  the 
size  of  CMM  for  these  curves,  the  width  of  the  pulse  used  in  the  moment 
method  region  was  kept  constant  as  0.05X. 


It  is  seen  from  the  results  in  Figures  3.8  to  3.10  that  as  the 
wedge  angle  or  the  wave  incident  angle  is  decreased,  the  dependencey 
of  the  surface  current  on  the  extent  of  the  moment  method  is  increased. 

Usually  the  6TD  solutions  for  the  surface  currents  on  the  wedge 
are  in  very  good  agreement  with  exact  solutions  except  for  the  region 
very  near  to  the  edge.  In  Figure  3.12,  the  GTD  and  our  hybrid  method 
solutions  for  the  currents  in  that  region  (i.e.,  up  to  0.05A  from  the 
edge)  are  compared  with  exact  solution  for  the  90°  wedge  and  incident 
angles.  It  is  seen  that  the  error  of  the  GTD  currents  becomes  large 
as  the  observation  point  approaches  the  edge  and  also  that  overall 
agreement  of  our  hybrid  method  solution  with  exact  solution  is  a  little 
better  than  that  of  the  GTD. 

All  the  results  in  this  chapter  were  obtained  with  the  second 
order  approximation.  Also,  by  experience,  the  second  order  approxi¬ 
mations  were  seen  to  be  enough  for  acceptable  results.  We  see  that 
fast  convergence  of  the  iterative  method  is  possible  with  the  appli¬ 
cation  of  this  method  to  a  planar  surface.  This  will  be  reinforced 
later  by  considering  a  square  cylinder  in  Chapter  IV. 
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♦  AS  -  MM 
XM  *  60  X 


d(X) 


Figure  3.5.  Currents  on  a  wedge  for  <j>  =90°, 
V30’  CMM=0'5x»  and  XM=60X.W 


■S 


34 


Figure  3.6. 
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♦  ♦♦  AS-MM 
X  M  *  50  X 


d(X) 

Figure  3.7.  Currents  on  a  wedge  for  <j>  *120°, 
V30  »  cmms°-5x*  and  *M=50X.W 


ra  3  8  Currents  on  a  wedqe  f or  y30 
re  _ toftO  r  * 0.5X,  and  XM-50X. 


AS  —  MM 


Figure  3.9.  Currents  on  a  wedge  for  4>  =30°, 
4>i*120°,  C^l.OX,  and  XM=50X.W 
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Figure  3.11.  Currents  on  a  wedge  for  $  =60°, 
<^=60°,  0,^=1. OX,  and  XM=40X.  w 
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CHAPTER  IV 


SQUARE  CYLINDER 


4.1  Integral  Equations 

The  geometry  of  a  perfectly  conducting  square  cylinder  in  two- 
dimensions  is  shown  in  Figure  4.1.  A  TE  plane  wave  is  incident  on 
the  square  cylinder  with  an  incident  angle  6^.  For  convenience,  we 

10  3  4 

represent  the  surfaces  AB,  BC,  CD,  and  DA  by  C  ,  C  ,  C  ,  and  C  ,  re¬ 
spectively.  The  moment  method  regions  are  set  up  around  the  four  edges 

12  3  4 

of  the  square  cylinder  and  are  denoted  by  C^,  CjjjM,  C^,  and  CjJjM, 

respectively  on  the  surfaces  c\  C2,  C2,  and  C4. 

It  is  assumed  in  this  geometry  that  the  surfaces  C1  and  C4  are 

2  3 

lit  and  that  the  surfaces  C  and  C  are  shadowed.  Therefore  the  in¬ 
cident  angle  ei  must  be  between  90°  and  180°  for  the  following  equa¬ 
tions  which  are  derived  below. 

The  currents  ^(R)  and  in  the  lit  regions  can  be  expres¬ 

sed  as 

^p(ff)  =  2n1xH.{ff)  +  2nTx/  ’  )xV 'G(r)dH 1  (4.1) 

C4 

and 
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(4.2) 


^p(tf)  -  ZfytfflO  +  2B4x  j  ^p(f)x7*G(r)d£'  . 

C1 

p  _  o 

Also  the  shadow  region  currents  (R)  and  (If)  can  be  written  as 

ur 


and 

^p(ff)  =  2*3xH\K)  +  2ft3x  I  f  ^p^'Jxv'GCrJdj,1  (4.4) 

k=l  p  k 

m  L 

where  n^,  n2,  n3,  and  n4  are  the  outward  unit  normal  vectors  to  the 
surfaces  C1,  C2,  C3,  and  C4,  respectively. 


The  first  order  approximate  currents  in  the  four  moment  method 
regions  can  be  obtained  from: 


3  3 

*  2n4x  I  f  l[jM(R,)xV,G(r)djl,+2n4x  l  f  ^p(f '  )x  VG(r)dt’ ; 
k=l  pk  k^2  pk 

V|HI  « 

MW 


Figure  4.1.  Surface  divisions  in  a  square  cylinder 
in  two  dimensions. 
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The  expressions  for  the  currents  in  the  asymptotic  regions  are 
the  same  as  Equations  (4.5)  to  (4.8)  except  that  the  observation  points 
are  moved  from  the  MM-regions  to  the  AS-regions. 

Observing  all  these  equations  we  can  see  that  the  observation 
and  source  points  are  not  on  the  same  surface  in  any  one  equation. 
Therefore  as  in  the  wedge  case  in  Chapter  III,  the  singularity  problem 
(i.e.,  exclusion  of  an  entire  patch  containing  the  singular  point) 
does  not  appear  in  the  square  cylinder  case. 

The  currents  J^R)  and  (E-)  can  be  solved  simultaneously  by 
the  iterative  method.  After  substituting  these  currents  into  the 

p  _  3  _ 

equations  for  ^p(R)  and  (R),  they  also  can  be  solved  by  the  same 
technique. 

4.2  Impedance  Matrix 

The  moment  method  is  used  to  solve  Equations  (4.5)  to  (4.8). 

For  this  purpose,  we  use  simple  pulse  basis  functions  in  the  point¬ 
matching  technique. 

As  we  expect  from  the  equations  we  derived,  the  resulting  im¬ 
pedance  matrix  will  be  very  complicated  because  there  are  four  unknown 
currents  involved  in  the  calculations.  Therefore,  we  will  obtain  an 
approximate  expression  for  the  matrix  by  ignoring  certain  terms  in 
the  equations.  For  the  sake  of  clarity,  expressing  Equations  (4.5) 
to  (4.8)  in  scalar  forms,  we  have 
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(4.9) 


The  result  of  the  differentiation  of  the  Green's  function  with 
the  associated  vector  calculation  is  expressed  as  which  is 

given  by 


VM')  -  -  4  H<2>(Br)  -i - 


(4.13) 
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If  we  substitute  Equations  (4.10)  to  (4.12)  into  Equation  (4.9) 
repeatedly,  we  obtain  an  equation  which  has  one  unknown  current  I^U) 
and  is  composed  of  multiple  integrals  containing  1^(1')  in  their 
integrands.  After  ignoring  the  triple  and  higher  multiple  integral 
terms,  we  have 

i  4  l 

IMMU)  =  4  k?2  \  ^  IMM^'  )Gq(®',,»^1  )Gq(®'»®'")d®’,dfc" 

CMM  CMM 


where 


+  Q(H)+2  J 


Q(**)GqU,r)dr 


(4.14) 


QU)  -  2  \  {  4pUMB,U,V)dr 

k-2  ck 


(4.15) 


Assume  that  the  solution  of  Equation  (4.14)  can  be  expressed 
in  the  series  form  of 


N 

\h 


P(*r*h) 


(4.16) 


where 


r 


1 


POD  =  < 


for  <-^ 


1^0  elsewhere. 


(4.17) 
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w 

and  m  =  -j^,  N  is  a  number  of  subintervals  in  W. 

Inserting  Equation  (4.16)  in  Equation  (4.14)  and  rearranging 
terms,  we  have 


=  QU)  +  2 


f  QU')G  U,i')dr 

c4 

lmm 


(4.18) 


Multiplying  both  sides  of  Equation  (4.18)  by  5(8,-^),  m=l,2,...N, 
integrating  them  over  the  surface  C^,  and  dividing  both  sides  by 
M,  we  have 


from  which  1^  can  be  determined  by  solving  a  matrix  equation  of  the 
type 
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f  Zmr^  ^nl  =  fVm^»  m,n=l,2,...,N.  The  elements  of  the  imped¬ 
ance  matrix  [Z^]  are  obtained  from 

4n  -  ’’(W  '  4  ],2/k  (4.20) 

CMM 

where  m,n=l ,2,3, . . . ,N. 

Also  the  elements  of  the  voltage  matrix  [V*]  are  determined  by 

Vm  =  Q(Am^  +  2  /  QU^O^,*')^'  m=l,2,...,N.  (4.21) 

r4 

In  Equation  (4.20)  the  most  dominant  term  is 

/  G  (n',an)G  Um,aMAJKfe' 

c4 

which  is  due  to  the  surface  current  on  the  lit  surface  cj^.  Therefore 
for  the  simpler  form  of  Z*n,  we  can  ignore  all  other  terms  except  the 
dominant  term. 
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4.3  Numerical  Results 


For  numerical  integration,  the  length  of  one  side  W  of  a  square 
cylinder  was  divided  into  60  subintervals.  The  moment  method  region 

L 

CjJjM(k=l,2,3,4)  was  set  up  with  one  half  of  W  in  it's  length.  Therefore 
the  extent  of  the  moment  method  region  (C^)  will  be  changed  according 
to  the  size  of  the  cylinder.  However,  as  we  have  seen  in  the  wedge 
problem,  the  current  on  the  surface  of  a  square  cylinder  also  will 
not  be  affected  by  C^M,  as  we  will  see  later  in  the  numerical  results. 

Numerical  calculations  were  performed  for  various  cylinder  sizes 
and  angles  of  the  incident  field,  and  comparisons  were  made  with  inde¬ 
pendent  GTD  solutions.  For  the  GTD  solutions,  up  to  triply  diffracted 
fields  were  considered,  and  the  diffraction  coefficient  formulas  ob¬ 
tained  by  Kouyoumjian  and  Pathak  were  used.  The  magnitude  and  phase 
of  the  surface  current  on  a  0.7X  square  cylinder  with  the  wave  incident 
angle  6^95°  are  shown  in  Figure  4.2.  It  is  seen  that  overall  agree¬ 
ment  with  the  GTD  solutions  is  very  good  except  that  the  phase  of  the 

3 

current  on  the  deep  shadowed  region  (C  )  deviates  slightly  from  the 
GTD  result.  Also  the  result  for  a  square  cylinder  with  W=1.2\and 
0^=120°,  which  is  shown  in  Figure  4.3,  is  seen  to  be  in  very  good 
agreement  with  the  GTD  solution. 

Results  for  larger  W,  that  is  W=3.0X  and  4.2X  are,  respectively 
shown  in  Figures  4.4  and  4.5.  We  see  that  overall  agreements  are  quite 
good  with  slight  deviations  in  the  deep  shadow  region.  Although  we 
obtained  results  for  square  cylinders  up  to  4.2\  on  a  side,  the  method 
does  not  show  any  limits  on  the  size  of  a  cylinder  or  on  the  wave 
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incident  angle.  Next  let  us  examine  the  results  of  the  surface  cur- 
on  a  3.0  square  cylinder  with  CMM=0.5A  and  2. OX,  respectively  in 
Figures  4.6  and  4.7.  The  number  of  pulse  basis  functions  in  is 

proportional  to  C^,  and  the  width  of  the  pulse  is  kept  constant  at 
0.05  for  those  data.  By  comparing  Figure  4.6  with  Figure  4.7,  we 
can  realize  that  the  larger  makes  the  result  slightly  better, 
especially  in  the  deep  shadow  region.  However,  it  is  clearly  seen 
that  the  extent  of  the  moment  method  region  does  not  affect  the  surface 
current  significantly. 

It  is  well  known  that  a  square  cylinder  can  have  an  interior 
resonance  problem  since  the  integral  equation  has  a  nonunique  solution 
at  the  resonant  frequencies.  The  moment  method  produces  an  erroneous 
solution  for  the  surface  current  on  a  square  cylinder  at  those  resonant 
frequencies.  The  geometries  of  the  square  cylinders  which  correspond 
to  the  resonant  frequencies  are  determined  by 


W  =  m2+n2  x 


m,n=0,l,2,3,. 


(4.22) 


from  which  we  can  find  an  Infinite  number  of  square  cylinders  which 
are  in  the  resonance  mode.  Some  of  these  geometries  are  obtained 
from  Equation  (4.22)  are  given  by  W=0.5x,  0.707\,  1 .OX,  1.118X,  1.5x, 

1.58X,  2. OX,  2.121a,  3.0a,  3.041A .  For  a  square  cylinder,  the 

resonant  frequencies  are  so  sharply  defined  that  we  may  not  detect 
the  resonance  phenomena  even  at  frequencies  deviated  very  slightly 
from  them. 
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Figure  4.8  shows  the  results  of  the  surface  currents  on  the 
0.707X  square  cylinder  with  0^=95°  which  is  one  of  the  geometries 
having  an  interior  resonance.  Three  methods,  the  GTD,  moment,  and 
hybrid  methods  are  used  to  obtain  those  results.  It  is  seen  that  the 
moment  method  produces  erroneous  results  in  both  the  magnitude  and 
phase  of  the  current.  Also  when  the  same  square  cylinder  with  0^=120° 
is  tried,  the  moment  method  solutions  become  worse  as  shown  in  Figure 
4.9. 

On  the  other  hand,  our  hybrid  method  solutions  do  not  show  any 
resonance  phenomena  and  instead  show  very  good  agreement  with  the  GTD 
solutions  in  both  cases  of  0^=95°  and  120°.  Results  for  W=1.118X, 
which  also  correspond  to  one  of  the  resonant  geometries,  are  given 
in  Figure  4.10.  Our  solution  does  not  show  any  resonance  effect. 

It  is  worthwhile  to  mention  that  the  GTD  and  our  hybrid  method 
solutions  for  W=0. 705  X (which  is  slightly  different  from  0.707X)  have 
remarkably  good  agreement  with  the  moment  method  solutions  as  shown 
in  Figure  4.11. 
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Figure  4.2.  Current  on  a  square  cylinder  for 
W=0.7X ,  0^95°,  and  C^O.SW. 
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PHASE  (DEGREES) 


0  1.2  2.4  3.6  4.8 

d(X) 


Figure  4.3.  Current  on  a  square  cylinder  for 
W=1.2X,  e.=120°,  and  CMM=0.5W. 
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Figure  4.6.  Current  on  a  square  cylinder  for 
W=3.0X,  ©1=1 10°,  and  0^=0. 5X. 
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4.7.  Current  on  a  square  cylinder  for 
W=3.0X,  0 .=1 10°,  and  CMM=2.0X. 


Figure  4.8.  Current  on  a  square  cylinder  for  W=0.707A 
(resonance  case),  0^95°,  and  CMM=0.5W. 


Figure  4.9.  Current  on  a  square  cylinder  for  W=0.707X 
(resonance  case),  0^=120°,  and  CMM=0.5W. 
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Figure  4.10.  Current  on  a  square  cylinder  for  W=1 . 1 18x 
(resonance  case),  e^llO0,  and  0^=0. 5W. 
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CHAPTER  V 
CIRCULAR  CYLINDER 


5.1  Integral  Equations 

We  will  consider  a  perfectly  conducting  circular  cylinder  which 
is  illuminated  by  a  TE  plane  wave  as  shown  in  Figure  5.1.  The  moment 

fl  c 

method  regions  which  are  represented  by  C^  and  C^,  respectively  in 
the  lit  and  shadowed  regions,  are  located  in  the  regions  near  the 
shadow  boundary. 


Figure  5.1.  Surface  divisions  in  a  circular  cylinder 
in  two  dimensions. 
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The  equations  for  the  surface  currents  on  a  circular  cylinder 
will  have  the  same  forms  as  the  general  equations  except  for  the  inte¬ 
gration  regions.  Therefore  we  can  write  the  following  equations  di¬ 
rectly  from  the  general  equations: 

^p(R)  =  2nxH1  (FT)  +  2^  x  /  J^p(R* )xVG(r)di'  (5.1) 

& 

Jop(R')  =  2nxH1(ff)  +  2n  x  /  J^p(R 1  )xV'G(r)d£' 

C* 


+  2n  x  /  ^qD(R  1  )xV'G(r)d£,' 

Cs  P 


(5.2) 


4<r) 


2n  x  /  lJM(R*)xV'G(r)dJi'+2fi  x  /  I„M(R'  )xv'G(r)d£' 

rs 

UMM  °MM 


+  2n  x  /  ^op(R ' ) XV ' G ( r ) cUt ' 
Cs 


(5.3) 


lmW 


=  2nx/  TL(R')*7*G(r)dr+2fl  x  j  I^jR'JxVGtrJdi 


'MM 


'MM 


(5.4) 


G(r)  is  defined  in  Equation  (2.6).  The  expressions  for  the  cur¬ 
rents  on  the  asymptotic  regions  have  the  same  forms  as  Equations  (5.3) 
and  (5.4)  except  the  observation  points  are  moved  from  the  moment 
method  regions  to  the  asymptotic  regions. 


64 


5 . 2  Consideration  of  the  Singularity  Patch 


As  we  mentioned  before,  the  equations  for  the  curved  surface 
have  as  many  singular  points  as  observation  points.  Therefore  Equa¬ 
tions  (5.1)  to  (5.4)  may  not  be  useful  for  practical  numerical  cal¬ 
culations  because  we  have  to  skip  some  finite  areas  containing  the 
singular  point  instead  of  one  point  for  the  evaluation  of  the  principle 
value  of  integration. 

Let  us  consider  the  exaggerated  part  of  the  curved  surface  shown 
in  Figure  5.2.  The  points  A,  B,  C  and  D  are  the  boundaries  of  di¬ 
vision,  and  P  ,  Pj  and  P2  represent  observation  points  which  are  lo¬ 
cated  at  the  center  of  each  division.  Thus,  the  currents  at  the  points 
Pp  PQ  and  P2  represent  the  currents  on  the  subintervals  AB,  BC  and 
CD,  respectively. 


0 


Figure  5.2.  Source  points  around  the  sir  ,.ar  point  PQ. 
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When  we  perform  numerical  integration  with  the  observation  point 


at  PQ,  the  contribution  of  the  source  currents  which  are  located  in 
the  interval  BC  will  be  disregarded,  and  thus,  the  errors  will  be  in¬ 
creased.  If  we  subdivide  the  interval  BC  into  many  subintervals,  then 
the  resulting  errors  will  be  decreased.  However,  this  way  of  solving 
the  problem  may  not  be  efficient  when  we  consider  limits  of  computer 
memory  storage  and  running  time.  Therefore  we  will  derive  the  magnetic 
field  integral  equation  which  will  be  useful  for  practical  calculations 
in  general  and  then  derive  the  specific  result  for  the  circular  cyl¬ 
inder  problem. 

Just  inside  the  surface  of  a  perfectly  conducting  body,  the  fol¬ 
lowing  boundary  condition  will  be  satisfied. 

nxTTi (R")  +  nxHs(W)  =  0  just  inside  S  (5.5) 

where  H1^)  is  the  incident  magnetic  field,  and  fls(ff)  is  the  magnetic 
field  due  to  the  electric  surface  current  on  S  and  is  given  by 


HS(R)  = 


V  x  /  J(ff')  G(r)  ds * 
S 


From  Equations  (5.5)  and  (5.6),  we  have 


(5.6) 


nxH*(ff)  +  nxyx  ^  J(ff' )G(r)ds '  =  0  .  (5.7) 

Using  the  vector  identity  7  x  (AB)  =  (7A)  x  B  +  A(7xl5)  and  noting  that 
7x3(1T')=0,  we  have 
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(5.8) 


r 


nxH1^)  +  tfxf  J(ff’)  x  7'G(r)ds'  =  0 
S 

where  we  used  the  relationship  7G(r)=-7'G(r) ,  and  V'G(r)  represents 
gradient  of  G(r)  in  the  source  coordinate  system.  If  we  divide  the 
surface  S  into  two  surfaces,  S-S^  and  which  is  finite  and  includes 
a  singular  point  at  the  center  of  it,  we  can  write  Equation  (5.8)  as 


nxHn(R)  +  nx 


/  J(R‘ )x7'G(r)ds*  +  nx  /  J(R' )xv'G(r)ds' 
S-Sf  S 


=  0. 


(5.9) 


Since  the  third  term  of  Equation  (5.9)  still  has  an  r"1  singu¬ 
larity  as  R  approaches  R',  we  will  consider  a  hemispherical  surface 
S£  which  is  mounted  over  the  singular  point  Pq  with  radius  rc  as  shown 
in  Figure  5.3.  Then  Equation  (5.9)  can  be  rewritten  as 


fixF^lT)  +  nx  /  3(ff')x7'G(r)ds’ 


S-S< 


+  lim  inx  /  U(R')x7'G(r)ds'  +  nx  /  J(fc' )xv'G(r)ds' 1=0 

V°l  VSe  Se  J 

(5.10) 


Let 


^■.(ff)  =  lim 

V° 


nx 


f  3(ff')x7'G(r)ds' 
S 

e 


(5.11) 
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Figure  5.3.  Surface  divisions  around  the  singular  point  P  . 


and 


J,(R)  *  lim  nx  f  )xV'G(r)ds'  . 

i  S  -0  s*-s 

e  f  e 


(5.12) 


Then  Equation  (5.10)  becomes 


ftxH1^)  +  fix  /  J(^')xV'G(r)ds'  +  "J^R)  +  ^(R)  =  0  .  (5.13) 


S-S 


Using  a  small  argument  approximation  for  H^(0r),  J^(R)  can  be  easily 
obtained  as  (see  Appendix  A) 


f 


Although  this  result  is  obtained  for  the  two  dimensional  problem,  the 
sarnie  result  can  be  obtained  for  the  three  dimensional  sphere  problem 
which  we  will  consider  in  Chapter  VI. 

Next  we  will  consider  Js(ft)  in  Equation  (5.12).  Let  be 
in  the  direction  of  £  on  the  surface,  then  we  can  write 

J(R')  =  t  J(W)  .  (5.15) 


The  unit  vector  t  will  be  -$  for  the  geometry  in  Figure  5.1.  Also 
let 


V'G(r)  =  r'  Gc(r) 


(5.16) 


where 


Gc(r)  s  "  J5  Hi2)(0r) 


(5.17) 


and  r'  =  — 5 - 

Since  f)  x  (t'x?')  =  t  cos0,  where  0  is  the  angle  between  the  unit 
vectors  n'  and  r'.  Equation  (5.12)  can  be  rewritten  as 


v 

?»■ 


71  r 
& 


£ 


k- 


X(ff)  *  lim  t  /  cos0  J(R’)Gr(r)ds’ 

£  se+o  sf-s  c 


(5.18) 
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or  since  n'*r'  =  cose,  we  can  write  Equation  (5.18)  as 

X(R)  =  1  im  €  f  n'-r*  J(R‘)GJr)ds'  (5.19) 

c  S  -0  S/-S  c 

e  t  e 

where  n'  is  the  unit  normal  vector  to  the  surface  at  the  source  point. 
To  obtain  the  approximation  for  ^(R),  we  assume  that 

J(R')  sj(lf)  (5.20) 

for  the  small  finite  area  Sf.  Then  Equation  (5.19)  becomes 

X,(ff)=3(ff)lim  f  Ti'-r'  Gr(r)ds*  .  (5.21) 

S>0  sj-s 

e  f  e 

From  Equations  (5.13),  (5.14)  and  (5.21)  we  have  expression  for  3(ff), 
which  is  given  by 


0(ff)  =  Cn  [2W(R)  +  2nx  f  J(E* )xV'G(r)ds'] 

S-Sf 


(5.22) 


where  the  constant  Cq  is 


Cn  = 


(5.23) 


1-2  lim  /  n'*r'  G.(r)ds' 

V°  srse 
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Now  let  us  calculate  Cq  explicitly  for  a  circular  cylinder.  Consider 
an  exaggerated  part  of  the  circular  surface  shown  in  Figure  5.4. 

Points  A  and  B  are  the  boundaries  of  the  finite  surface  S^,  and  the 
observation  point  P  is  located  at  the  center  of  AB.  The  value  of  n'*r' 
can  be  obtained  from  Figure  5.4  as 

n'*r'  =  cos  0  =  Yg  .  (5.24) 


Figure  5.4.  Surface  around  the  singular  point  P  to  give 
relationship  between  a,  r,  and  6 . 


Also  the  small  argument  approximation  for  Gc(r)  is  given  by 


71 


Sc(r)  =  -  |j  (  +  j  ~)  for  small  gr. 


(5.25) 


Substituting  Equations  (5.24)  and  (5.25)  into  Equation  (5.23),  we 
obtain 

Cq  * - fr - .  (5.26) 

1-4  y+o i  (^)("  + j  ^}dr 

r(fu  r0 

Evaluating  the  integration  in  Equation  (5.26),  we  obtain 


Notice  that  for  very  small  rc>  Cq  becomes  one,  and  thus.  Equation 
(5.22)  becomes  the  usual  magnetic  field  integral  equation. 

The  Equation  (5.22)  will  be  useful  for  practical  numerical  inte¬ 
gration  because  the  finite  area  is  removed  instead  of  just  a  point 
for  evaluation  of  the  principle  value  of  integration.  The  constant 
Cq  will  supplement  the  contribution  which  is  lost  by  removing  the 
finite  area  Sf.  Although  the  result  in  Equation  (5.27)  is  obtained 
for  a  circular  cylinder,  the  same  result  may  be  used  as  an  approxi¬ 
mation  for  an  arbitrary  smooth  surface  in  two  dimensions. 
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5.3  Impedance  Matrix 

To  solve  for  Imm(R)  a"d  T*,(R).  we  use  the  moment  method.  For 
the  sake  of  simplicity,  we  will  derive  the  approximate  form  of  the 
impedance  matrix  by  ignoring  all  possible  terms  that  make  a  minor  con¬ 
tribution.  Expressing  Equations  (5.3)  and  (5.4)  in  the  scalar  forms, 
we  have 


!&,(*)  ■  2  /  l5H(R')Gd(ia')dr  -2  f  lSM(lt')Sd(R.ir')d*' 


r  4  t  n 


S  /TT  i 


rH 


'MM 


and 


ImM(R)  =  2  /  IMM(*')Gd(R’*,)dr+2  *  I5M(R,)Gd^')dr  (5*29^ 


S  /77 1 


'MM 


'MM 


where 


G„(iuf')  ■  -*fLSzrLHi2)(8|R-R,|) 


(5.30) 


Substitution  of  Equation  (5.29)  in  Equation  (5.28)  leads  to 


ifiMw  *  ^  f  4<I'>Gd(R>f')dr 


r$. 

°MM 
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4  /  /  4(R,)Gd(R,,^,)6d(R^,,)dii,dr 


"'MM 


CMM 


(5.31) 


If  we  substitute  Equation  (5.29)  into  Equation  (5.31)  repeatedly,  we 
will  have  an  infinite  number  of  multiple  integral  terms.  For  an  ap¬ 
proximation,  we  ignore  all  multiple  integrals  higher  than  double  in¬ 
tegrals.  This  approximation  can  be  justified  by  noting  that  the 
dominant  terms  in  Equation  (5.31)  are  the  first  and  last  terms. 
Therefore  Equation  (5.31)  can  be  reduced  to 
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(5.33) 


Let  the  current  I$M(R)  be  approximated  by 


where 


r 


P(R)  = 


1 


for  |R-R„|  <  £ 


0  elsewhere 


(5.34) 


where  AS.  is  the  length  of  one  subinterval.  Substituting  Equation 
(5.34)  into  Equation  (5.32)  and  rearranging  terms,  we  have 


N 

Ji 


p(R-ffn)-2Gd(R\Rn)A£  -4  f  Gd(R '  ,*n)Gd(R,ff'  )AJLdS  1 


"MM 


*2f  Jop(R,)Gd(R.R,)d4,  •  (5.35) 

CS 


Multiplying  both  sides  of  Equation  (5.35)  by  <S(R-tfm),  m=l,2,...,N, 
.  .  0 

integrating  them  over  CMM,  and  dividing  them  by  M,  we  have 


i/. 


P(VRn)  -  2  MW** 


-4/ 


Cs 

UMM 


Gd(R'  ,Rm)®d(Rm,R' 
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(5.36) 


If  we  represent  Equation  (5.36)  in  the  form  of 

KinlllJ]  »  [Vfjl  m,n,=l,2,...,N,  (5.37) 

we  have  the  impedance  matrix  [Z^]  whose  elements  are  given  by 

ZL  *  P<VV  -  2  GA-V“ 

-  4  /  Gd(R'^n)Gd(Rm,R')Aild£'  .  (5.38) 

Cs 

lmm 

Also  the  elements  of  (V^l  are  determined  from 

*i*  2  5  4(f')Sd(V)<lf  •  (5-39) 

Cs 

5.4  Numerical  Results 

We  divide  the  circumference  of  the  cylinder  into  120  sub¬ 
intervals  for  the  various  numerical  calculations  which  follow.  The 

D  Q 

total  extent  of  the  MM-region  on  the  lit  surface,  C^,  was  60  in  angle, 
and  20  pulse  basis  functions  in  the  point-matching  technique  were  used 
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K 


. - 

to  obtain  the  current  in  the  MM-region.  Our  results  (marked  by  +  in 
all  data)  were  compared  with  the  exact  so  utions  (marked  by  the  solid 
lines).  Usually  the  second  or  the  third  order  approximations  were 
enough  for  acceptable  results.  Also  all  results  were  obtained  with 
the  constant  Cq  given  in  Equation  (5.27),  which  was  noted  by  CqJ*1  in 
the  figures  to  follow. 

Figures  5.5  to  5.8  show  the  magnitude  and  phase  of  the  surface 
currents  on  the  circular  cylinders  of  0.2X  to  5.7X  radii.  It  is  seen 
that  agreement  with  the  exact  currents  are  quite  good  regardless  of 
the  cylinder  size.  By  comparing  Figure  5.9  with  Figure  5.10,  which 
show  the  current  results  on  a  5.2X  circular  cylinder,  respectively 
with  C^jM=30°  and  cj^=120°,  we  can  see  the  effect  of  the  extent  of  the 
MM-region  on  the  surface  currents.  It  is  shown  that  the  surface  cur¬ 
rents  are  affected  very  little  by  the  size  of  C^jM,  which  is  a  very 

desirable  feature  of  this  method.  However,  it  will  be  seen  that  as 
£ 

Cmm  is  increased,  a  slightly  better  result  is  obtained. 

The  effect  of  the  constant  CQ  is  shown  by  Figures  5.11  to  5.14. 
The  results  in  Figures  5.11  and  5.12  were  obtained  with  a  circular 
cylinder  of  the  4.8X  radius.  We  see  that  a  little  improvement  was 
made  with  the  constant  Cq.  However,  for  the  radius  of  3. OX,  the  ninth 
order  approximation  gave  the  results  in  Figures  5.13  and  5.14,  where 
it  is  clearly  seen  that  the  constant  Cq  provides  a  great  improvement 
in  the  magnitude  and  a  little  improvement  in  the  phase  of  the  surface 
;  -  current.  We  realize  that  the  radius  of  4.8X  corresponds  to  a  non- 

r 

$  resonant  frequency  and  that  the  radius  of  3. OX  is  near  a  resonance 

i 

f 

* 

% 
f 


n 


mode.  For  the  geometries  corresponding  to  the  resonant  frequencies, 
the  higher  order  approximations  are  shown  to  be  needed  to  obtain 
acceptable  results. 

Recently,  several  methods  such  as  the  combined  source,  combined 
field,  and  hybrid  methods  have  been  presented  to  solve  the  interior 
resonance  problem.  However,  it  seems  that  any  of  these  methods  are 
difficult  to  apply  to  our  method  because  we  use  the  iteration  technique 
for  7*  (10  and  3*p{R),  which  needs  a  proper  initial  starting  current. 

Actually  the  interior  resonance  does  not  affect  the  surface 
current  on  a  circular  cylinder  so  seriously.  Although  we  obtained 
the  results  for  the  cylinders  of  up  to  the  radius  of  5.7X,  the  method 
does  not  have  difficulties  for  the  larger  cylinders.  Also  there  are 
no  difficulties  in  this  method  for  the  cylinders  of  radii  smaller  than 
0.2X. 
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Figure  5.8.  Current  on  a  circular  cylinder  for 
a=5.A  ,  CjL=60°,  and  C 
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Figure  5.11. 


Current  on  a  circular  cylinder  for 
a=4.8  X»  C|^=60°,  and  CQ=1. 
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9 ( DEGREES ) 


Figure  5.12.  Current  on  a  circular  cylinder  for 
a=4.8X,  cL=60°,  and  CJl. 
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Figure  5.14.  Current  on  a  circular  cylinder  for 
a=3.0X,  0^=60°,  and  C^l. 
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CHAPTER  VI 
SPHERE 


6.1  Integral  Equations 

So  far,  we  have  considered  conductors  only  in  two-dimensions 
where  integrals  over  the  surface  reduced  to  line  integrals.  In  this 
chapter,  we  will  consider  a  three  dimensional  problem,  that  of  a  per¬ 
fectly  conducting  sphere. 

Since  the  integration  must  be  performed  over  the  surface  of  the 
sphere,  a  longer  computer  running  time  and  increased  memory  storage 
will  be  needed  for  the  three-dimensional  problem  than  for  the  two- 
dimensional  one.  However,  by  using  the  symmetric  property  which  a 
body  of  revolution  has,  we  can  decrease  the  computing  time  very  greatly 
and  also  the  computer  memory  storage  needed. 

The  general  equations  we  derived  in  Chapter  II  cannot  be  reduced 
to  a  simpler  form  for  a  sphere,  but  instead  it  will  be  useful  to  derive 
integral  equations  for  the  surface  currents  in  the  E-  and  H-planes. 

The  geometry  in  which  a  sphere  is  illuminated  by  a  plane  wave  is  shown 
in  Figure  6.1.  Let  the  incident  magnetic  field  be  xpolarized  and 
traveling  to  the  negative  z  direction  as  shown  in  Figure  6.1.  At  an 
observation  point  P  the  induced  surface  current  'J(R’)  can  be  separated 
into  two  components  QJq(^)  and  $J^(1T),.  where  0  and  £  represent  the 
unit  vectors  in  the  directions  of  0  and  $,  respectively. 
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Figure  6.1  Geometry  of  a  sphere  illuminated  by  the  x-polarized 

plane  wave. 

Consider  the  magnetic  field  integral  equation  given  by 


J(ff)  •  ^(K)  +  2nx  f  W)xV'G(r)ds' 
S 


(6.1) 


where 


TJ1  (K)  =  2nxH^(R)  . 


(6.2) 


Since  ^(R)  and  U(TT)  are  composed  of  the  §  and  $  components,  they  can 
be  written  as 
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(6.3) 


n 


^00  =  0  J^ff)  +  $  J^ff) 

and 

700  -  $J0OO  +  fyOO  .  (6.4) 


The  currents  J0OO  and  J(j)(TT)  can  be  obtained,  respectively  by  9*700 
and  $•‘300,  which  are  given  by 


Je(R)  «  J0OO  +  20«n  X  f  J(T?‘  )xV'G(r)ds' 

(6.5) 

J<j,(R)  =  J4OO  +  2$*n  x  £7(10 )xV'G(r)ds' 

(6.6) 

A  A  A 

By  using  the  vector  identity  A-IxIT  =  and  noting  that  6xn=-$  and 

$xn=@  on  the  spherical  surface,  we  can  rewrite  Equations  (6.5)  and 

(6.6)  as 

O0OO  =  Jg(T0  -  2  /$ »7(T0' )xV'G(r)ds' 

(6.7) 

J.(10  -  jJ(10  +  2  /  §*7(10 )xV'6(r)ds' 

♦  ♦  5 

(6.8) 

Substituting  Equation  (6.4)  into  Equations  (6.7)  and  (6.8), 

we  obtain 

JQ00  =  jj(10  -2  /  $-0'J0(R')44'J^')1xVG(r)ds' 

S  J 

(6.9) 

f 

. 

i 

:  ! 

\ 
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v*}  a  J«}W  +  2  {  [s*e,J0(^H4,J4(^,)]^,6(r)ds'  («*io) 

where  §'  and  $'  are  the  unit  vectors,  respectively  in  the  0  and  <t> 
directions  at  a  source  point  P'.  Let 


V* G(r)  *  r'6s{r) 


(6.11) 


Gs(r)  -  -  ( J0  +  7)  -4irr 


(6.12) 


Using  the  notations  given  by 


CAfl(r)  =  $•§'  *  r'  Mr) 


(6.13) 


CflA(r)  =  $.$'  x  G Jr) 


(6.14) 


Vr)  =  x  Gs(r) 


C..(r)  =  @4'  *  Gs(r) 


(6.15) 


(6.16) 


we  can  express  Equations  (6.9)  and  (6.10)  as 


y*>  *  -  z{  [wh<*'>  +  Vr,vHds'  (6J7) 


J4(D-jJm*2f[Vr)J#(f)+.CH(r)J^)]df  .  (6.18) 
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After  calculating  the  vectors  In  Equations  (6.13)  to  (6.16),  we  can 
obtain  the  following  results 


t 


C00(r)  =  C^(r)  s7  ^ine,sine+cos(4,-4,')(cosecos6'-l)jGs(r) 

(6.19) 
and 

C0^(r)  =  -  C  (r)  =  f  (cose-cose  ')sin(4>-$')Gs(r). 
where  r  is  given  by 

r  =  a  J2^1-sinesin0,cos(^-(|),)-cosecosel| 
and  a  is  the  sphere  radius. 

_ -j  ^ 

Since  H  (ff)  is  x-polarized,  J0(lf)  and  J^(T0  are  given  by 

J^(ff)  =  2  sirxfr  eJ8acos8  (6.22) 

J^(1T)  =  2  cose  cos4>  e#^050  .  (6.23) 

Equations  (6.17)  to  (6.20)  are  not  very  useful  for  practical 
numerical  calculations  since  so  many  observation  points  are  involved 
in  them.  Therefore  we  will  derive  Integral  equations  having  only  the 
E-  and  H-plane  currents  as  unknowns.  For  this  purpose,  if  we  represent 
the  E-  and  H-plane  currents  of  ^(1T)  as  «^(e)  and  0^(e),  respectively. 
Equations  (6.22)  and  (6.23)  can  be  expressed  as 
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(6.20) 


(6.21) 


.1 


1 


jJ(ff)  -  Jj(0)sin4>  (6,24) 

jJ<*)  =  Jj>(0)c°s«j)  (6.25) 

where  J0(e)  and  J^( e)  are  independent  on  the  angle  4  and  are  given 
by 


Jq(  0)  =  2  ej **cos 0  (6.26) 

J^(0)  =  2  cos 0  eJ Sacos 0  (6.27) 

Also  the  relationships  of  Equations  (6.24)  and  (6.25)  are  valid  for 
Jiff)  and  Jiff).  Thus, 

J0(ff)  *  JQ(0)  sin^  (6.28) 


J^(ff)  -  J^e)  cos*  .  (6.29) 

Equations  (6.28)  and  (6.29)  are  useful  properties  which  can  be  applied 
to  an  arbitrary  conducting  body  of  revolution  when  IT^ff)  1$  x-polarlzetj. 
Using  Equations  (6.24),  (6.25),  (6.28)  and  (6.29)  in  Equations  (6.17) 
and  (6.18),  we  can  obtain 

JQ(0)  *  Jj(0)  -  2  /  [c1(r1)J0(01)  +  C^r^J^'jjds* 
and 

V0)  *  j4{0)  +  2  {  [Dl(r2)Je(0,)  +  ^2 ^ r2 ) ^ ® ' )] ’ 


(6.30) 

(6.31) 


I 

i 
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where 


w  • 

c2(rl)  - 
°l(r2)  = 

= 


a®s ( r i )  r  1 

— — —  ^s1n0sin0'+s1nV  (cos9cos9'-l)jsin<|>' 


aGs(rl)  2 

- prp  ( cos  9-cos  e* )  cos' <t>* 

aGs(r2)  ? 

- r —  (cos 6-cos e1  )sin  <j>‘ 


aGs(r2)  f  1 

- r —  |^inesine,+cos<$>'(cosecose'-l)jcos<l), 


(6.32) 


(6.33) 

(6.34) 


(6.35) 


Also  and  r2  are  given  by 

r-j  =  a  J2  (l-sinesine’sin^'-cosecose')  (6.36) 

r2  =  a  J2  (l-sinesine'ccs^’-cosecose1)  .  (6.37) 

If  the  currents  Jo(0)  and  J.(e)  are  known,  the  surface  currents  at 

0  <P 

any  arbitrary  point  on  the  sphere  can  be  readily  determined  from  Equa¬ 
tions  (6.28)  and  (6.29).  Equations  (6.30)  and  (6.31)  can  be  applied 
directly  to  any  equations  we  derived  in  Chapter  II  by  adding  proper 
subscript  and  superscript. 

If  we  look  at  Equations  (6.32)  to  (6.37),  we  can  easily  see  that 
D](r2)  and  D2(r2)  can  be  obtained,  respectively  from  C2(rj)  and  C^r^) 
just  replacing  0'  by  90°  -  4)'.  This  characteristic  will  make  the  com¬ 
puter  program  simple. 
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6.2  Consideration  of  the  Singularity  Patch 

For  reasons  we  discussed  in  Section  5.2  of  Chapter  V,  we  need 
to  consider  the  singularity  problem.  We  will  use  a  part  of  the  result 
of  Section  5.2  to  derive  the  proper  constants  for  the  sphere.  If  we 
rewrite  Equation  (5.13)  here  again  with  the  different  definition  of 
V'G(r)  given  by  Equation  (6.11),  the  equation  will  be 

tfxH1^)  +  n  x  /  3(ff' )xy'G(r)ds'  +  ^(R)  +  ^(R)  =  0  (6.38) 

S-Sf 

where  J^E)  and  ^(ff)  are  given  in  Equations  (5.11)  and  (5.12),  re¬ 
spectively. 

By  the  similar  procedure  as  done  in  Section  5.2,  we  can  easily 
find  T| (R)  in  the  three  dimensional  case,  which  is  given  by  (see  Ap¬ 
pendix  B) 

(ft)  =  -  ^  <J(R)  .  (6.39) 


Now  let  us  derive  ^(ff)  for  the  @  and  $  components.  Using  the  results 
of  Equations  (6.30)  and  (6.31),  we  can  write  J20(e)  and  as 

J2e(e)  -  -  2  Hm  J  jc^r^fo')  +  C^r^  (e  1  )J  ds*  (6.40) 
eT  f  ^  g 

and 


V9) 


2^0  s/s 

e>  f  e 


(6.41) 
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where  ^(e)  and  ^(0)  are  the  currents  observed  on  the  <J>=90°  and 
4>s0°  lines,  respectively. 

Since  the  surface  area  is  small,  we  can  make  the  following 
approximations  for  J2Q(e)  and  <^(9) 

<t>‘  ~  7  on  Sf  for  J20(0)  (6.42) 

<J>'  ~  0  on  Sf  for  ^(e)  •  (6.43) 

Then  using  Equations  (6.42)  and  (6.43),  the  approximations  for  C^r^), 
C2(>"])»  D-|(r2),  and  D2(r2)  on  Sf  can  be  obtained  as 


CT(rT)  a  -  2!  Gs(ra)  on  Sf 

(6.44) 

^2(^1 )  ~  0  on  Sf 

(6.45) 

and 


D^)  =0  on  Sf 

(6.46) 

r. 

°2(r2) =  ’  2a  Gs(raJ  on  Sf * 

(6.47) 

where 


(6.48) 


ra  =  ay2(1-sinesine '-cosecose ' )  on 

Substitutions  of  Equations  (6.44)  and  (6.45)  in  Equation  (6.40)  and 
Equations  (6.46)  and  (6.47)  in  Equation  (6.41)  yield 

J2  (6)  =  -  2  lim  }  (-  ^f)Gs{ra)J  (e‘)ds-  (6.49) 

o  -+U 

e  '  e 

J2*(0)  “  2  lira  /  (-^.(rJJJe^ds'  .  (6.50) 

S  +0  s;-s  s  a  $ 

e  f  e 

Also  for  the  small  surface  area  S^,  we  can  make  the  following  approxi¬ 
mations. 


J0(e')  -  J0(e)  on  sf 

(6.51) 

J  (e‘)  =  J.(e)  on  sf  . 

<p  1 

(6.52) 

Then  by  substituting  Equations  (6.51)  and  (6.52)  into 

Equations  (6.49) 

and  (6.50),  we  obtain 

-r 

o2e(0)  =•-  2  J6(e)  H<nQ  J  (-55)Ss(ra)ds' 

e  f  e 

(6.53) 
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(6.54) 


X 


9 

i 


V9’  « 


J.  (© ) 

<P 


Tim 
S  +0 


S#-S 


(■?f>Gs<ra>ds' 


Assuming  that  the  patch  is  a  circular  surface  with  the  radius  rg, 
then  we  can  obtain  the  following  results  from  Equations  (6.53)  and 
(6.54). 


033(6)  a-  SQ  0Q(e)  (6.55) 

^2c^0^  ~S0  ^  (6.56) 

where 

S0  -?m[2  -  <Z^Jers)e'Jers]  .  (6.57) 


Then  from  Equations  (6.38),  (5.11),  (5.12),  (6.55)  and  (6.56)  we  can 
obtain  the  following  integral  equations. 

1  (6.58) 

I .  -  TT 
J4>  j 


1  (6.59) 

J<()  *0 


JL(0)  =  Sfl  0  •  |  2HXH1  (ff)+2nx  /  3(ff' )xs"G(r)ds 

0  0  L  s-sf 


0(e)  =  $  *  \  2nxH1(ff)+2?ix  /  3(f')x7'G(r)ds 

♦  *  l  S-Sf 


or  by  using  the  results  of  Equations  (6.30)  and  (6.31),  we  can  express 
Equations  (6.58)  and  (6.59)  as 


V0)  *  V 


Jq(0)  -2  $/s  [c1(r1)Je(0-)+C2(r1)J4)(0')jds|  (6.60) 


V0)  *  s<j, 


{J>> t2  4 


Dl(r2)Je(e*)+D2(r2)J4,(e 


*  )j  Js^  (6. 


61) 


where 


c  -  1 

0  _TFSS 


(6.62) 


S<fr  =  *  (6*63) 

We  see  that  as  r.  approaches  zero,  SQ  and  S,.  approach  one.  This  means 

b  y  <p 

that  if  the  sphere  surface  is  divided  into  very  many  subpatches  then 
the  constants  and  will  not  affect  the  integral  equations  sig¬ 
nificantly.  However,  dividing  the  surface  into  many  divisions  may 
not  be  efficient  for  numerical  calculations.  By  using  the  constants 
^  and  we  can  provide  the  same  effect  as  using  a  large  number  of 
divisions  for  numerical  Integration. 

6.3  Impedance  Matrix 

In  calculating  the  Impedance  matrix,  the  moment  method  region 
is  set  up  around  the  shadow  boundary  making  a  ring  form  as  shown  in 
Figure  6.2.  If  we  consider  all  points  on  as  unknowns,  the  computer 


100 


storage  may  not  handle  the  matrix  because  of  too  many  unknowns.  Ther 
fore  we  will  use  the  results  we  have  derived  in  Section  6.2  to  obtain 
the  impedance  matrix. 


Figure  6.2.  Surface  divisions  in  a  sphere. 


Also  since  the  equations  for  currents  are  involved  with  two 
currents  J0(e)  and  J^(0),  the  expressions  for  the  impedance  matrices 
Zeinn  and  Z  will  be  very  complicated  and  may  not  be  useful  for  nu¬ 
merical  calculations.  Therefore  we  will  derive  comparatively  simple 
approximate  forms  by  ignoring  all  the  terms  that  contribute  insignif- 


=  ®  *  W*> 


(6.64) 


and 


(6.65) 


Then  we  can  write  integral  equations  for  IMe(0)  and  1^(0)  directly 

from  Equations  (6.30)  and  (6.31).  The  currents  1^(6)  and  I^(e)  rep¬ 
resent  the  E-plane  and  H-plane  currents,  respectively  in  the  moment 
method  region.  For  the  sake  of  clarity,  we  will  use  the  notations 
C1(e;0‘,<J>'),  C2(O;0 * ,<|> * ) ,  D1  (e ;e * ,<f> ’ )  and  D^eie',*')  instead  of 

cl(r]).  Cgfrj),  O^rg),  and  02(r2),  respectively.  Then  ljg(0), 
ijSl  ^0 ) ,  Ij^O).  and  1^(0)  can  be  expressed  as 


+  |(0)  (6.66) 
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where 


/  [d^  «  e' ,  V  )iM0(e'  )+D2(e;e 
JMM 


+  2 


/  [D1(e;el^,)i5e(0,)+D2(9;e,*<t>,)IM^0,)]ds’ 

SMM 


+ 

<P 


(6.67) 


)+c2(0;e'»4> 


ds* 


-  2 


](0»9'  »4»'  )^mq(0)+C2(0»0'  »4>i  )I 


^(9)  3  2 


+  2 


[D1(0;0,,<D,)l5e(01)+D2(0;0l^')I^(e') 


MM 
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and 


F0(e)  =  -  2  finite;  €f,e *)j^ (e,)+c2(e;e',4>l)j^)(e,)J ds1  (6.70) 


F»(e)  =2  f [°i (e ;e ' 4 1  )v3^q (e ' )+D2(e ;e  1 )  J^(0 ' )1 ds 1  (6.71) 

cS  J 


JP9 (e ' )  and  Jp$(0')  rePresent  the  9  and  ♦  components  of  Jgp(e'), 
respectively.  Substituting  Equations  (6.68)  and  (6.69)  into  Equations 
(6.66)  and  (6.67)  repeatedly,  and  ignoring  the  third  and  the  higher 
multiple  integrals,  we  obtain 

IMe(e)  =-2/A[ci(0;e,.<t>,)IMe(e,)+c2(0;e,.*,)IM4>(0, )]  ds' 


Hr! 


+4  /  f  I 

SMM  SMM 


Me 


(0*)  C1(0";0’ ,<t>’ )C1(e;e",q>")-D1  (0" ;0 ' ,4> '  )C2(0; 0" ,<|>" )J  ds'ds" 


4/$  /  ^  )^i (0*6M»<|>M) “D2 ( 0 " ; 0 ‘  a q> '  )C2(0» 9" ><!>")  ds'ds" 


s3  ^ 

^MM  bMM 


+  F0{®) 


(6.72) 


tfyobZf''  [D1(0;0,„')IMV0')+D2(0;0,,,.)IMy0')]ds' 


Hr! 


%  \  40(e,)[■C^(0,l;0,’♦,)D^(9;e,,^,,)+Dl(e,,;0,^,)D2(0;0,,*♦M)Jds,ds,, 


S  S 

Wyy  Jyy 

HH  Hr! 
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H  f  f  IM(J)(e,)[-C2(e";0,'♦,)D1^0i0^<l>")+O2(9“50^<t•,)D2(eiO^4'")jds,dS,, 
Syy 

nn  Pin 

+  l  (e)  (6.73) 

Substituting  Equation  (6.73)  into  Equation  (6.72)  repeatedly, 
and  then  ignoring  the  triple  or  the  higher  multiple  integrals  contain¬ 
ing  I*(e')  and  the  double  or  the  higher  multiple  integrals  containing 
F^Ce),  we  obtain  the  following  approximation  for  IjSj0(e). 

^e*0)  *  -  2  /  IM0(e’ )C-|(e;e* »<*>' )ds' 

^MM 

-4/  f  ij5j0(0')Di(eH;e,,(|>,)c2(e;0,',<|,")ds'ds" 

SMM  SMM 

+  4f  /  XM0^®'  )Ci(0"l0'  »♦'  )C1(9;0",(|,M)ds,ds" 

+  Fe(0)  -  2  /  F^(0')C2(e;0,,^,)ds’  .  (6.74) 

si, 

riFi 

Also  by  a  similar  procedure  we  can  obtain  an  expression  for 
!^(e)  as 
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IW6>  -  2  f  lM,><e'>D2(e;e\4,')<is' 


+  F4>(0)  +  2/  Fe(e,)O1(0;e*.<t>*)ds'. 

4 

Assume  that  lj^(0)  and  lj^(0)  can  be  approximated  by 


I*-  / 


where 


P(0)  -  < 


II 

p<e-«n> 

N  . 

„I, 

■Jr  p(«-8n> 

r 

|  1 

for  |0-0n|<  f- 

0 

elsewhere. 

(6.75) 

(6.76) 

(6.77) 

(6.78) 


and  A0 


=  ^  ,  where  N  is  a  number  of  divisions. 


Substituting  Equations  (6.76)  and  (6.77)  into  Equation  (6.74),  and 
rearranging  terms,  we  have 


L 


<P(e-e„) 


+ 


2/  c1(e;en,*,)<r 


+  4/  f  D-j  (6*1  ;©n,<|)  *  )C2(6;0M  ,<|>")  d<t>  *  ds " 
SMM  ^MM 


-4/  /  C1(0";en,<|.,)C1(0;0",4>,,)d<(.,ds"  1  a2sinenA0  ' 

SMM  *MM 


Fg(0) 


(6.79) 


where  the  integration  region  <j>^M  represents  the  lit  moment  method 
region  for  the  angle 

Multiplying  both  sides  of  Equation  (6.79)  by  6(0-6m),  m=l,2,...,N, 
integrating  them  over  0^,  and  then  dividing  both  sides  by  A0,  we 
obtain 
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+  4  f  f  D1(eI1;0n,4.')C2(em;0,l,*")d<J.1ds" 

SMM  ♦mm 

-  4  i  f  C1(0^0n,4>,)C1(0m;0%r)<rds"]  a2sine^e  r 

^  MM 
nn  nn 


W-2l 

4 


(6.80) 


Representing  Equation  (6.80)  in  the  matrix  form  of 


"•n-1-2 . N 


(6.81) 


means  that  the  elements  of  the  impedance  matrix  [Z^]  are  given  by 


Z*  = 
©mn 


p(V6n>  +Y  5 

^  *MM 


+  4  f  /  d1  (e";0n»4>’ )c2(em;e",<t>")d  'ds" 


SMM  ^MM 


-  4  /  f  C1(e'';en,^')C1(0m;0",<j.,')d4,,ds" 


Ss  J- 


a2sin0nA0  .  (6.82) 
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Also  the  elements  of  [V*^]  can  be  obtained  from 


(6.83) 


By  the  same  procedure  we  can  obtain  the  impedance  matrix  and 
the  voltage  matrix  for  the  $  component  of  the  surface  current,  whose 
elements  are  given  by 

2U  -  p«Ven>  +{-  2  i  D2<Ven-*'>d*' 

^  ^MM 


+  4f  f  c2(e- ;en.* •  )D, (©nse- )d4» ' ds- 


(b  ini 

MM 


-4f  f  D2(e";en,(i.,)D2(9m;9,,,<)>,,)ct|>,ds"|'a2sin0nA0 
9s.  ) 


5mm  *mm 


(6.84) 


and 


'4m  *  F*<9m>  +  2  f  Fe Ce •  )Di (em:e ■  •  )ds ■  .  (6.85) 

SMH 

The  dominant  terms  in  Equations  (6.82)  and  (6.84)  are  the  first 
two  ones.  Therefore  if  we  want  the  simpler  forms  for  computer  pro¬ 
grams,  we  may  ignore  the  last  two  terms  in  those  equations. 


6.4  Numerical  Results 

Two  possible  segmentation  schemes  for  numerical  integration  over 
the  sphere  surface  are  shown  in  Figures  6.3  and  6.4.  In  the  uniform 
segmentation  shown  in  Figure  6.3,  the  surface  is  divided  into  equal 
angles  in  the  8  and  $  directions.  Therefore  the  patch  areas  around 
8=0°  will  be  very  small  compared  to  those  around  8*90°. 


X 


Figure  6.3.  Uniform  segmentation  scheme. 


On  the  other  hand,  in  the  variable  $  segmentation  shown  in  Figure 
6.4,  we  divide  the  surface  equally  in  the  9  direction  and  proportional 
to  sinQ  in  the  $  direction  so  that  all  patch  areas  are  almost  equal. 
Experience  shows  that  the  better  numerical  results  are  obtained  with 
the  variable  $  segmentation  scheme.  All  numerical  results  for  a  sphere 
were  obtained  with  the  variable  q>  segmentation  scheme,  where  we  divided 


1 


•4 


X 


Figure  6.4.  Variable  $  segmentation  scheme. 

the  surface  of  the  interval  6=0  to  n/2  into  twenty  four  bands  (we 
call  the  number  of  this  division  N)  and  each  band  starting  at  the  pole 
(9=0)  in  the  interval  <JfO  to  tt/2  was  divided  into  subpatches  as  fol¬ 
lows;  1,2,4,5,7,8,10,11,13,14,15,16,18,19,20,20,21,22,22,23,23,24, 
24,24.  With  the  larger  number  of  N,  better  results  can  be  obtained, 
but  also  the  resulting  computational  time  will  be  increased  signifi¬ 
cantly.  For  example,  the  computer*  running  time  with  N=24  is  about 
twenty  six  minutes  for  the  seventh  order  approximation,  but  if  we 
increase  the  number  N  by  two  times,  the  computing  time  will  be 
increased  by  approximately  six  times  for  the  same  order  approximation. 


♦Modified  Oatacraft  model  6024. 
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This  example  shows  that  there  is  the  upper  limit  on  the  sphere  size 
which  can  be  handled  practically  by  this  method. 

The  interior  resonance  appears  to  affect  the  surface  current 
more  seriously  on  a  sphere  than  on  a  circular  cylinder.  This  seems 
to  be  due  to  the  interaction  of  the  two  currents  J0(0)  and  J^(e). 

Also  it  is  observed  that  the  resonance  problem  becomes  more  severe 
for  the  larger  spheres.  Therefore  this  resonance  phenomenon  limits 
the  numerical  calculations  for  arbitrary  sizes  of  a  sphere  but  may 
not  be  of  any  consequence  for  other  3-dimensional  geometries. 

We  obtained  acceptable  numerical  results  for  conducting  spheres 
of  radii  0.2&  to  2.3A  using  N=24.  Compared  to  a  circular  cylinder, 
higher  order  approximation  were  needed  for  a  sphere.  The  surface 
current  on  a  sphere  was  obtained  with  the  seventh  to  the  eleventh  order 
approximations.  Also  in  the  moment  method  region,  whose  size  was  30° 
in  angle  for  all  data,  eight  pulse  basis  and  delta  weighting  functions 
were  used.  All  results  were  compared  with  the  exact  eigenfunction 
solutions  which  were  denoted  by  solid  lines  in  the  figures. 

In  Figures  6.5  and  6.6  the  H-plane  and  E-plane  currents,  which 
are  denoted  by  and  JQ,  respectively,  are  shown  for  a  sphere  of 
radius  a=0.25A.  The  agreements  with  the  exact  solutions  are  remarkably 
good  in  both  the  E-plane  and  H-plane  currents.  Also  the  results  for 
a=0.5A  are  in  good  agreement  with  the  exact  solutions  as  shown  in 
Figures  6.7  and  6.8. 
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Figures  6.9  and  6.10  show  the  surface  currents  for  the  1.0X 
sphere,  where  we  can  see  that  the  differences  from  the  exact  solutions 
are  larger  in  both  the  magnitude  and  phase  for  the  surface  current 
compared  to  the  results  obtained  for  the  smaller  spheres. 

It  seems  that  the  surface  current  on  a  sphere  is  affected  not 
only  at  the  resonance  frequency  but  also  the  frequencies  around  it. 
Therefore  we  cannot  find  many  sizes  of  spheres  with  which  we  can  obtain 
numerical  results  by  this  method.  For  example,  although  a=1.5A  and 
2. Ox  are  not  the  exact  geometries  that  have  the  interior  resonance, 
the  surface  currents  on  those  spheres  are  much  different  from  the  exact 
currents.  The  radii  of  the  spheres  which  are  in  the  interior  resonance 
modes  around  a=1.5X  and  2. OX  are  a  1.48A  and  1.98x,  respectively. 

The  surface  current  results  for  the  sphere  of  1.7x  are  shown  in  Figures 
6.11  and  6.12  and  are  in  fairly  good  agreement  with  the  exact  results. 
However,  it  is  obvious  that  they  are  not  that  good  compared  to  the 
results  for  the  spheres  of  0.5A  or  1 .OX  radius. 

The  largest  radius  of  a  sphere  for  which  we  obtained  numerical 
results  was  2.3X,  and  the  results  are  shown  in  Figures  6.13  and  6.14. 

We  can  observe  that  even  though  the  results  are  acceptable,  the  mag¬ 
nitude  curves  show  many  oscillations  compared  to  other  sizes  of  spheres. 
The  greatest  departure  from  the  exact  solutions  is  shown  in  the  phase 
of  the  H-plane  current  in  the  shadow  region.  The  oscillations  shown 
on  the  magnitude  curves  can  be  decreased  a  little,  but  not  signifi¬ 
cantly  by  increasing  the  number  of  divisions.  However,  as  it  has 
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been  mentioned  before,  by  increasing  the  number  of  divisions  the  com¬ 
puting  time  will  be  greatly  increased.  Therefore  instead  of  increasing 
N,  it  will  be  more  efficient  to  use  the  constants  ^  and  which  are 
given  in  Equations  (6.62)  and  (6.63)  to  make  improvement  in  the  re¬ 
sults.  Notice  that  the  results  in  Figures  6.13  and  6.14  were  obtained 
with  those  constants. 

If  we  examine  the  results  shown  in  Figures  6.15  and  6.16,  where 

we  used  N=24  without  the  constants  SQ  and  S. ,  we  can  realize  that 

u  <p 

those  constants,  which  were  derived  on  the  basis  of  very  careful 
approximations,  make  the  magnitude  curves  considerably  smoother. 

Although  we  obtained  numerical  results  for  spheres  of  up  to  the 
2.3X  radius,  it  will  be  possible  to  obtain  the  surface  currents  for 
larger  spheres  by  using  a  larger  number  of  divisions  with  the  constants 
Sq  and  . 


114 


0 


45.0 


90.0 

6  (DEGREES) 


135.0 


180.0 


0 


-  EXACT 

+  ++  AS-MM 


Figure  6.5.  The  H-plane  current  on  a  SDhere  for 
a=0.25X,  N=24,  and  SMM=30  . 
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Figure  6.6.  The  E-plane  current  on  a  sphere  for 
a=0.25A,  N=24,  and  5^=30°. 
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Figure  6.7.  The  H-plane  current  on  a  sphere  for 
a=0.5x,  N=24,  and  SMM=30°. 
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Figure  6.8.  The  E-plane  current  on  a  sphere  for 
a=0.5X,  N*24,  and  5^=30°. 
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Figure  6.10.  The  E-plane  current  on  a  sphere  for 
a=l .OX,  N=24,  and  ^MM=30°. 
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Figure  6.11.  The  H-plane  current  on  a  sphere  for 
a=1.7X,  N=24,  and  S^-300. 
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Figure  6.12.  The  E-plane  currents  on  a  sphere  for 
a*l.A,  N=24,  and  S^^O0. 
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Figure  6.13.  The  H-plane  current  on  a  sphere  for 
a=2.3x,  N=24,  S^O,  and  S0t  S^l. 
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Figure  6.14.  The  E-plane  current  on  a  sphere  for 
a=2.3X,  N=24,  5^=30°,  and  Sq.S^I. 
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Figure  6.15.  The  H-plane  current  on  a  sphere  for 
a=2.3A,  N-24,  S^O0,  and  S0tS^«l. 
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Figure  6.16.  The  E-plane  current  on  a  sphere  for 
a=2.3 N*24,  5^=30°,  and  SQ,S<|)=1. 


CHAPTER  VII 


SUMMARY  AND  CONCLUSIONS 


We  have  derived  the  necessary  integral  equations  for  the  hybrid 
AS-MM  method  and  applied  them  to  four  different  shaped  perfectly  con¬ 
ducting  bodies,  namely  a  wedge,  square  and  circular  cylinders  and  a 
sphere. 

The  modified  geometrical  optics  current  terms  which  were  denoted 

_ 0  _  _ c  _ 

by  Jgp(R)  and  JQp(R),  respectively  in  the  lit  and  shadow  regions  were 
obtained  by  an  iterative  technique  in  numerical  calculations.  The 
currents  in  the  moment  method  regions,  which  were  denoted  by  Tj^(R) 

C  M 

and  ( R ) *  resPective^y  in  the  lit  and  shadow  regions  were  solved 

by  a  matrix  technique  ( i . e . ,  the  moment  method).  For  this  purpose, 
we  derived  the  impedance  matrix  for  each  conducting  body.  The  currents 
in  other  than  the  moment  method  region  were  determined  by  using  the 
currents  in  the  MM-region  in  the  magnetic  field  integral  equation  for 
the  first  order  approximation.  For  the  higher  order  approximation, 
the  iterative  method  was  used  to  obtain  both  the  moment  method  and 
asymptotic  region  currents.  The  extent  of  the  moment  method  region 
did  not  affect  the  surface  current  significantly.  Usually  the  second 
or  the  third  order  approximation  produced  good  results  for  the  wedge, 
square  and  circular  cylinders,  but  for  the  sphere  higher  than  the  fifth 
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order  approximations  were  required  to  give  acceptable  results.  This 
is  maybe  due  to  the  interior  resonance  problem. 

The  singularity  patch  problem  in  the  curved  surface  geometries 
was  considered  and  some  constants  were  derived  on  the  basis  of  approxi¬ 
mations  to  make  up  some  of  the  contribution  lost  by  excluding  the  patch 
area  instead  of  just  the  singular  point  in  the  calculations.  It  was 
seen  that  the  constants  made  possible  good  improvement  in  the  results 
of  the  surface  current  and  that  they  were  useful  for  practical  numeri¬ 
cal  calculations. 

In  this  study  we  considered  a  technique  for  combining  a  moment 
method  current  and  an  asymptotic  current  in  both  two  and  three  dimen¬ 
sions.  From  the  formulations  for  the  surface  current  on  a  conducting 
body,  we  see  that  the  method  does  not  need  any  a  priori  knowledge  of 
the  current  form  away  from  the  moment  method  region.  Also  we  don't 
have  to  follow  ray  paths  on  the  surface  of  a  conductor  as  we  do  in 
the  GTD  method.  Therefore  these  advantages  of  this  method  tend  to 
make  it  relatively  easy  to  obtain  the  surface  current  on  arbitrarily 
shaped  bodies. 

As  we  have  seen  from  the  several  examples,  this  method  does  not 
have  difficulties  in  its  application  to  the  arbitrary  sizes  cf  bodies. 
However,  this  method  usually  requires  relatively  large  computing  time 
and  memory  storage  especially  for  bodies  having  many  surfaces.  This 
disadvantage  may  limit  the  method  in  practical  numerical  calculations. 


As  we  have  seen  in  the  wedge  problem  in  Chaper  III,  this  method 


cannot  provide  us  acceptable  results  for  the  surface  current  on  the 
wedge  if  the  magnetic  incident  field  is  grazing  to  the  lit  surface 
or  if  the  wedge  angle  is  small.  Of  course  these  limitations  are  of 
no  consequence  in  bodies  whose  dimensions  are  finite. 

Also  it  was  shown  that  this  method  could  not  eliminate  what  ap¬ 
pears  to  be  spurious  resonances  appearing  in  the  large  sizes  of  a 
sphere.  However,  this  method  was  not  affected  by  the  resonance  in 
its  applications  to  the  square  and  circular  cylinders. 

In  spite  of  these  disadvantages  and  limitations,  this  method 
provides  generally  good  results  for  the  surface  currents  on  the  per¬ 
fectly  conducting  bodies.  Extension  of  this  method  to  more  complex 
bodies  should  be  possible.  Future  work  should  study  ways  to  reduce 
the  computer  running  time  by  making  appropriate  approximations  in  the 
integrals.  Future  work  should  also  consider  the  possibility  of  com¬ 
bining  this  hybrid  method  with  the  GTD.  Thus,  one  would  treat  areas 
far  removed  from  shadow  boundaries  and/or  sharp  discontinuities  in 
geometry  by  the  GTD  thereby  eliminating  numerical  integration  time. 
The  hybrid  method  of  this  paper  could  then  be  used  to  handle  a  wide 
variety  of  discontinuities  in  geometry  for  which  diffraction  coeffi¬ 
cients  are  unknown. 
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Appendix  A 


CURRENT  CALCULATION  AT  A  SINGULAR  POINT  FOR 
TWO-DIMENSIONAL  PROBLEMS 


To  derive  Equation  (5.14)  from  Equation  (5.11),  we  consider  a 
semicircular  path  i as  shown  in  Figure  Al.  Let  the  Surface  current 
■J(IT)  be  in  the  direction  of  £,  then  Equation  (5.11)  can  be  expressed 
as 


TUff)  =  t  lim  /  J(ff')  (n'-r1)  d*  (Al) 

e+O  "e 

where 

G(e)  =  -|HJ2)(8e)  .  (A2) 

For  small  0c,  H^(ge)  can  be  approximated  as 


H<2)(0e)=  1  -  j  \  loge  (Xjt) 


(A3) 


where  y  is  Euler's  constant. 


Then 

3e 


can  be  evaluated  as 


36(e)  _ 
3e 


1 

"  £rre 


(A4) 


o 


Figure  .  Al.  Integration  path  around  the  singular  point 
P  in  two  dimensional  problems. 

Also  for  small  circular  path  S,£,  J(E')  can  be  approximated  as 
Substituting  Equation  (A4)  into  Equation  (Al),  and  noting  that  n'-r'=l 
and  J(ff'  )=J(TT) ,  we  obtain 

=  £  J(l£)  lim  /  ( 

€*0  1 

€ 

~  J(J£)  lim/  (-  d0 

s+0  o 

a  -  7  J(ff)  (A5) 

which  leads  to  Equation  (5-14). 
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APPENDIX  B 


CURRENT  CALCULATION  AT  A  SINGULAR  POINT  FOR 
THREE-DIMENSIONAL  PROBLEMS 

Consider  a  hemispherical  surface  Se  with  the  radius  shown  in 
Figure  B1.  By  the  same  way  as  done  for  Equation  (Al),  Equation  (5.11) 
can  be  written  as 


J.(R)  =  t  11m  Jj(R‘)(n’-r*)  ds'  .  (Bl) 

1  0  S£ 

Using  -  -  (j8  +  -^)  ,  ds  '*  £?sin  0'd 0* d 4>* ,  and  approximation 

J(IT)  for  small  surface  ^  in  Equation  (Bl),  we  have 


2iro/2  •»  _  -  j$c  9 

J.(f)  =  t  J(E)  lim  J  J  (-jB-  e  sine'de'di))' 

1  e+0oo 


7(ff)  lim  2Yj2  (j%~)  ej0e  sine'de'd*' 
e-*0  oo  f 


2mr/2  • 

W) 

0  0 


W) 


(B2) 


which  leads  to  Equation  (6.39). 
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